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PREFACE. 


An abstract of the following investigation was publish ed in Nos. 
(2024-2026) of tbe Astronomisahe Naohriclden. But, a few verbal in- 
accuracies, together with numerous typographical errors, being found in 
the work as printed, it was decided to reprint the whole investigation, 
giving the mathematical developments more in detail, and also develop- 
ing the formula for the latitude of tlie moon, a subject not touched 
upon in that investigation. The author hopes to find leisure for the 
complete development of the perturbations of the moon’s motions, by 
means of the differential equations given in this first chapter. It is 
believed that these equations will give the co-ordinates of the moon in 
a more direct and simple manner than any combination of the general 
differential equations of motion hitherto employed for that purpose. 


iii 




Theory of the Moon’s Motion. 


INTBODUCTION. 

It is now nearly two centuries since the theory of the universal gravitation 
of matter according to the law of the inverse square of the distances was dis- 
covered and subjected to calculation by Newton. The theory when regarded as 
a law of nature may be stated as follows: Every particU of 'matter m the U7iivcrse 
attracts every other partide rvith a force which varies directly as its masSj and 
imerTidy as the square of the distance between them, 

Matlicmaticians have attempted to deduce from tins law of matter all the 
phenomena attending the motions of the heavenly bodies. For this purpose the 
Bun, planets and satellites may be regarded as particles of matter when compared 
with the distances which separaki them from the other bodies of the system, the 
masses of the planets being regarded as infinitely small in comparison with that 
of the sun. The principal mathematical consequences resulting from the opera- 
tion of tins law of matter may bo stated as in the three followii\g theorems : 

L The orhits of the planets and eo^mets arc conic sections m which the sun 
ocevpics the pyincipal focus; 

IL The radius vector of each pkmd or comet sweeps over equal areas in equal 
times; and 

III. The squares of the tmes of revohitmv of the different qda^icts are to each 
other in the same proportion as the cubes of their mean distances from the sun. 

These three laws, wliich were disciovcred by Keplee, may bo regardixl as the 
embodiment of the theory of universal gravitation, and are the foundation of 
physiiral or mathematical astronomy. They are, however, exact only on the sup- 
position that the masses of the planets are infinitely small in comparison with 
that of the sun, and hence arc not npplioablo to the solar system witbout some 
modification, because tbo masses of some of the planets are finito instead of being 
infinitely small. The motions of the planets thex’eforo do not strictly conform to 
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the preceding laws , for, in obeying then mutinl atti action, tlioy must (h m lU a 
little from the elliptical paths which they woul<l t\ icily follow il they wcu 
attracted only by the sun Not onlv ne tlio pi mots distuibul in tluii motions 
around the sun by reason of their mutual aiti action, but the satellites lu iRo 
distuibed in their motions around then piimiiies by the itti leiion ol the sun nul 
by the other planets of the system The sm illness ol tlu i^liiu t ii> in i s,s nx 
comparison with that of the sun, howcvci, pi units tin sc l)c:)(lu s to eon 1mm \ciy 
nearly to the laws of the elliptical motion, and hence tlu t ihul ition of th ii ju i 
tuibations is not a difficult pioblcm But in the motions of tlu s iP llib s ibmit 
then primaries, the principal distuibing body is the sun , anil tin c (li 1 1 ol lus 
attiaction is too great to be overlooked Espeeiilly is tins the ( i-i m n to 
the motion of tho moon aiounil the eaith, the di^tuibiiii^ fmu ot tlu un l)i ing 
so great as to cause the ti ansveisc axis of he i oibit to tie si i ibe a comj>li ti n vo 
lution in about nine years, mel also causing tho notles ot the oibit lei xcVoUc in 
about double that inteival 

2 The mathematical theory of the pcituibitions of tlu pi nets and ittlhlos 
develops the laws according to which tho vnums (c^rces opei ile in tlu di tmiuiiei 
of each other s movements They m ly be cl issod is follows 

First The opeiation of a comp n itively luge luiee duiin^ i '-hoit intiii il 
of time, 

Second The opeiation of a comp natively snull fotu dining i Ion ^ pi nod 
of time, and 

Third A very slow change in the mean motions of tin dilb u nt 1 ndn ^ ui ing 
fiom the variations of tho elements ot then mbits piudund b) tin u mulinl 
attractions 

The first class pi oduccs the ponodic incfpulitns, whnh d» pnul f lil\ on iIm 
mutual distances and coiifigui itmns of tlm dilU i< nt bruin tin *-1 * rnnl f 1 1 pio- 
duces the inequalities of long pci lorl, wlm h mi} 1 m s ml tr> (1» |m nd un tin runh^u- 
lation of the elements of thou motions, ind tlu lliiid <1 1 piodin tin r ^ ul u 
inequalities The second and thud dissrs iffrct onK tin rm ui mull 

3 The problem of dodueing ill tin ciKumlum wlmb if] 111 n U* n ut 

the different bodies of the system, due r tly fiom tin pnmiplt f»l ni u i i1 i t\i 
tation, has, poihips, moie thin any ollit i td tin in nin \ m n i 1* nu 

ticians and astronomeis dming the 1 ist c( ntiu v ind i In 1 1ln t!i i < <1 i!ic 

periodic and secular inoqiulities of tho j>l uir t \\« it dt v* lujad m tb ♦ i t ti t ul 

for the purposes of astronomy about the close of tlu li I ctntinj Con nh i ibh 

progress was also made m the deaclopmc nt of tlio Inn ii tlu t.m b\ f \ T r m i , but 
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tlie magnitude of tlie disturbing forces being so much greater than those which 
disturb the planetary motions, the labor necessary for the construction of a 
perfect lunar theory is vastly increased. La Place states that ho has determined 
all the inotpialities of the first, second and third orders, and the most impoitant 
ones of the fourth order, and has continued the approximation to (piantities of 
the fourth order inclusively, and also retained those of the fifth order, which arise 
in the calculation. But much more elaboi’ate theories of the moon s motion were 
published during the fu’st half of the present century, in which the approxima- 
tions are carried to terms of a much higher order than was attempted l>y 
La Vlxce] the most complete and perfect, being that of Plana, in which the 
approximations are carried to terms of the seventh order, was published in the 
year 18S2. Pontecoulant also investigated the lunar theory in the fourth 
volume of his Tlieoric Analiiti(][UQ du Sy^ithne du MondCj which was published in 
the year 184G, in which he has also carried the approximations to terms of the 
seventh order. The theories of the moon’s motion by Plana and PoNTKConLANT 
wore productH.l by entirely difibrent methods of developmeiit—— the one by employ- 
ing tlie true longitude as the imlejieudcnt variable, and the other using the mean 
longiimle for tlic same purpose — and when reduced to the same lorni were found 
to bo almost idcuiically the same, the coefficients of the iueipialitics seldom dif- 
fering from cae.h other by so much as one second ol arc. It would theroforo 
Bocm that the caleulaticms of l)oth had been correctly made. But porhiips tlio 
most elabornite ami conudeto investigations on the tlicory of the moons motion are 
those of Hanskn nml Delaunay. Hansen pulfiished very eluhorato do 

la Lam" in the ycfir 1857; but the aiiulvsis by whicb lie obtained his ibrraulio 
has not, so i;ir as I am informed, yet been published. The lunar theory by 
Delaunay was pulilished in two large rjnarto volumes, which appeared in the 
years 1850 and 1857, ami is ptuliaps tlie most complete and perfect worlc on the 
lunar tln^ory iha,t has yet been givmi to astronomers. 

4 * In a valuable and interesting discussion on the apparent inequalities of long 
period in tlu’ moiion of the moon, ]Tul>lishe<l in the Joiivnol of Srumoc 

and A rh (or September, 1870, Prof ‘ssor Newcomb makes the following statement 
in regard to the imxjiinlities of sliorl. period in tlie motion of the moon : 

*‘Tlie prolilern of (h4<u’mining ilie motion of the moon nr(')nnd the earth undoi 
the influence c>f the combined nilraelinn of the sun and planets has, more than 
any oilier, called forth the efforts of matliernaticians and astronomers. Lr^aily 
every grr‘at geometer since Nfavton has added something to the sirnplic^ity oi the 
acoAiracy of the solution, and in our own day wc luive seen it successfully cour- 
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pleted in its simplest form, in which the eaith, the moon and the sun are regarded 
as mateiial points, moving under the influence of their mutual attractions The 
sati'^factory solutions are due to the genius of Hansen and of Delaunay Work 
mg independently of each other, each using a method of his own invention moie 
rigorous than had before been applied, they arrived at expressions for the longi 
tude of the moon which, being compared, were found to exhibit an aveiage dis 
crepancy o± less than a second of arc Ho doubt could remain of the substantial 
correctness of each ’ 

This statement, coming from so eminent an astronomer, would lead us to 
suppose that the mathematical theory of the moon s motion, in so far as it depends 
on the attraction of the sun and planets, may be regarded as perfect Four years 
later the same writer, after a thorough and extended examination of the question, 
assures us of the correctness of his conclusions 

Hansen has also given in his tables of the moon two inequalities of long 
period depending on the action of Venus Delaunay, in computing the coef 
ficients of these two inequalities by his own method, reproduces, very nearly, 
Hansen s coefficient of one of them, but finds the other to be insensible 

Professor Newcomb then proceeds to inquire whether we have in either theory 
a satisfactory agreement with observations, and finally arrives at the conclusion 
that the problem of the inequalities of long period in the moon s mean motion is 
really no nearer such a solution as will agree with obseivation than when it was 
left by La Place , and he is obliged to resort to the operation of irregular and 
extraneous causes, the data for the accurate determination of which are wanting, 
in order to fill up the gap between theory and observation 

In the Monthly Notices of the Jtoyed AstTonoTniccd Society for November 14, 
1873, Mr Airy also enters into a discussion in regard to the inequalities of long 
period m the moon’s mean motion, and states his conclusion in these words I 
conceive it to be totally impossible that a complete and correct theory can leave 
such large and systematic discordances And I express my opinion that there is 
some serious defect in the Lunar Theory ’ 

These statements, coming from such eminent astronomers, are sufficient to 
encourage the belief that the published theoiies of the inequalities of long period 
in the moon s mean motion possess but very little value 

6 Let us now inquire whether the inequalities of shoit period in the moons 
motion are any more perfectly represented by existing theories than the inequali- 
ties of long period The theories of Hansen and Delaunay, according to Pro- 
fessor Newcomb, may be regarded as identical in their results, and they also 
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contain all tlie important researches of astronomers up to a very recent date. 
The theories of Plana and Pont:ecoulant may also be regarded as identical in 
their results; the former was published in 1832, and "TaSZes of the Moon," 
founded on it, were published by Professor Peiece in 1853. Plana’s theory was 
somewhat modified by the theoretical investigations of Hansen ; and some 
empirical corrections, proposed by Aiey and Longsteeth, were also introduced ; 
so that the tables in their present shape do not exhibit the results of pure theory 
alone. But if the empirical corrections introduced are legitimate parts of a cor- 
rect lunar theory, and only await the demonstrations of the physical astronomer, 
we may assume that they have a permanent value, and that the tables are perma- 
nently improved by their introduction ; otherwise they possess only a temporary 
value, and will vitiate the tables alter* a while to the extent of their primitive 
improvement. Now, the ephemeris of the moon, contained in the American 
JEphemeris and Mautioal Ahnanac, is derived from Peirce’s tables of the moon, 
and wo have only to consult the volumes of the Naval Observatory at ‘Washington 
in order to learn to what extent tho moon deviates from the requirements of 
theory. In this comparison we may suppose that the moon’s mean place is cor- 
rectly given by tho tables, in whicdi ca.se the observed deviations between theory 
and oliservation will servo to indicate the magnitudo of one or more inequalities 
of short period in the moon’s motion which analysis has yet failed to point out. 
Wo shall tlicroforo find that, in tho year 1856, the moon’s place differed from tho 
computed place between tlie limits + 10".4 and - 12".4; in 1857 tho errors were 
between tho limits +7".7 and -10".4; in 1858' they were between H14".4 and 
-9",3; in 1861 they were + 10".6 and “9".6; in 1862 they wore +9".0 and 
— 11".3; in 1871. they wore -18".4 and -~4".8; and in 1872 they were H- 14",0 
and — 12".0. These residuals of coin-se show only the oltscrved deviation of tho 
moon from her computed place, and do not niujessarily indicate tho maximum 
extent to which the theory deviates from the true place of tho moon, as cloudy 
weather no doubt frequently prevents observations at times when tho deviations 
are at their maximum values. Tho residuals above given would seem to indicate 
that tho mean place of the moon was very closely given by tho tables, since the 
positive and negative errors are very necwly equal to each other. According to 
Professor Newcomb, Hansen’s tables represented the moon’s place at tho time of 
their puVdication and for a few years after much more closely than Peiece’s. In 
tho year 1862 the errors of Hansen’s tables wome included within tho limits of 
+ 4".0 and — 9".0; in 1871 the limits were + 2".5 and -10".4; while in 1872 
they wore ™ 0".2 and — 16".2 ; thus indicating a narrower range of periodic 
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variations, and also a much, wider range in a progressive mean motion It should 
he obseived, however, that the preceding comparison exhibits only eirors of light 
ascension, and it is probable that the errors in longitude would generally be some 
what larger 

If we in like manner compare the observed and computed longitudes of the 
moon as given in the Seductions of the Greenimeh Lunar Observations fwm 1831 
to 1851, after correcting for the error of mean longitude, we shall hnd an oscilla- 
tion between theory and observation amounting to an average value of about 16" 

La Place says, in the Introduction to his Tlieoyy of the Moon, that the erior 
of the tables formed from his theory will very rarely evceed 32" , and he also 
states that the astronomer Bueg, by deriving the forms all the arguments fiom 
theory, and rectifying the coefficients bv means of numerous observations, had con- 
structed tables of the moons motion whose greatest errors were less thin 13" 

If this statement of La Place, which was made three quarters of a century ago, 
was borne out by the observations of the moon made at that time, it would 
indicate a degree of perfection in the lunar tables, with lospect to tho inequali- 
ties of short period, which is scarcely exceeded by the tables in use at the present 
tune 

6 Prom the preceding comparison it would appear that tho lunar theories of 
Hansen and Delaunat, when measured by the ciitcnon to which all physical 
theories must be subjected, are but little, if any, moie perfect than the theories 
of Plana and Pont]5cottlant, which preceded them by twenty years , and we 
also perceive that the lunar theory employed in tho Seductions of the Qieenwich 
Lunar Observations, and also the tables in use at the close of the last cc ntury, in 
so far as the inequalities of short period are concerned, are but little inferior to 
the tables in use at the present time It is also possible that tho supmoiity 
of recent tables over those in use half a century ago may ansc more fiom the 
corrections to the mean elements of her motion than from any nnpiovimcnt in 
the theory of her perturbations It would thus appear eithci th it tho theory of 
gravitation had been imperfectly or incorrectly developed, or else th it the moon 
undergoes perturbations from the action of other forces than gravitation 

7 The residuals above given would seem to indic ate th it our present lunar 
tables, instead of being correct to terms of the seventh order, nro really erroneous 
by some of the smaller terms of the third order, and as the writer had, piovious 
to making this comparison, assured himself, by a careful oxaminition of the 
mathematical theory of her motion, that some terms of tho thud order h id been 
overlooked, he does not hesitate to announce the fact to astronomers, and he 
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confidently believes that a correct theory developed to terms of the fifth order 
will be found to represent the motions of the moon, in so far as the inequalities of 
short period are concerned, with far greater precision than any published at the 
present day. 

8. It is true that wo have the combined assurance of all the great mathema- 
ticians of the present century who have given especial attention to the subject 
that no inequalities have been overlooked, and that the theory possesses all the 
accuracy claimed for it. To this we would reply that the error is fundamental 
and })rocedos any development of the perturbing function, simply growing out of 
the latitude of the perigee in the development of the undisturbed elliptical 
motion. If this statement is correct, it is easy to perceive that the fundamental 
error would iieeossarily ramify the whole development of the perturbing function 
and vitiate more or less all the conclusions deduced from it. Nor do we regard 
the argument that the agreement of so many profound mathematicians, that the 
theory of gravitation when legitimately applied to the moon’s motion would pro- 
duce the existing lunar theories, as possessing much weight. Such evidence is 
nmrely negative in its chfiractcr, and possesses no importance whatever when 
coming in conflict with positive evidence to the contrary. If any number of per- 
sons develop any mathematical expression by as many different methods, and all 
obtain identically tlie same result, it is satisfactory proof that the development 
has been (X)rrocUy made ; but it is no moro certain or satisfixetory than it would 
be if a single person had made all the developments by as many dillcrent methods. 
Besirlos, if an equation is supposed to possess a certain physicial or geometrical 
property, the development of the equation into an infinite series, or in any other 
form whatever, must also possess the same property. But no correct development 
or transformation of an equation can introduce any now element, either physical 
or geometrical. The development must reproduce all tho errors or Infirmities of 
the function from whicli it was derived. The assumption that a given mathe- 
matical expression possesses a certain physical or geometrical property should 
therefore b{^ fully justified by a rigorous discussion before such expression is made 
the basis of a physicial theory. It is to this point in the lunar theory tliat we 
would now call attention, 

IK For tliis purpose let us take tho fundamental equations for the latitude 
and rcraprocal of the projected radius vec.!tor of tho moon which are uhcmI by 
La and Plana as the basis of their respc^ctivo theories of. the moon’s 

motion. If wo denote tho mean distance and eccentricity of tho moon’s orbit by 
a and e, the longitudes of tho perigee and node by <o and Q, tlio inclination of the 



Xll 


THEORY OP THE MOONS MOTION 


orbit to the fixed plane by and the longitude, latitxide and radius vector by 
Vj 0 and r, and also put 

y = tan s ~ tan (1) 

we shall have the following equations, in which u denotes the reciprocal of the 
projected radius vector 

5 = y sin('y — Q), (2) 

«=— I— fl/r^+ecos(v-fi^)}=-i- (3) 

a(l — r) TQOSd 

La Place and Plana have supposed that equation (3) is the equition of a pro- 
jected ellipse, and we shall now proceed to show that it is not such an equation 
except for the particular case in which the transverse axis lies in the plane of 
projection 

Equation (2) gives 

l/l + s^==l/l +/sm^(v--Q) =i/l + tan^^ = — ^ (4) 

cos 0 

Whence, cos 8 = ^ , (5) 

l/l + y^ siTo? (v “ Q>) 

Bind: 

l/l + y* Bin* (i> — Si) 

Now, since u= — I—, we stall evidently have the maximum and minimum 
rcos 0 

values of u when r and cos 6 are respectively a minimum and maximum Now, 
the minimum value of r in an ellipse is r = a(l — e), and equation (5) gives 
1 

coad = — when cos^ is a minimum, and the maximum values of r and 

V'l+y* 

eoad are a(l+e) and 1, respectively, we shall therefore find 



maximum value of , 

a(l-e) ’ 

( 7 ) 

and 

minimum value of u= 

a(l H e) 

(8) 


If we now substitute the value of \/l + s* in equation (3), it will become 


— ^ ^ - 1 7-^ -IT I t/ 1 + y* sin*(u — Q) + e cob(u — a») } , 

rcoad a(l — e*) '■ ' 

and this equation gives 

maximum value ot u = ^ , 

a(l-e*) ’ 


( 9 ) 

( 10 ) 
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and minimum value of u = . (11) 

It will thus be seen that while the two minimum values of u agree, the two 
maximum values differ ty ’the quantity -(iey®) very nearly, or by a quantity of 

the third order. We shall soon see that the reason of this coincidence of the 
minimum values arises from the fact that the transverse axis of the orbit is in 
the plane of projection, in which case equation (9) is correct. 

If we now multiply equation (9) by cos 9 = we shall find 

■j/ 1 + 7 %in'\'y — Q) 

1= ^/l + e< iOBiv-w) ) 

r a(l — e®) [ -j/i -|- ain\v —Q)} 

In an ellipse the maximum and minimum values of r are a(l + e) and a(l— c), 
respectively. If then in equation (12) we suppose that to = we shall find the 
maximum and minimum values of r to be a(l + e) and a(l-e), as they evidently 
should be ; but if we suppose the ellipse to revolve in its own plane through an 
angle of 90°, <a will become SI + 90°, and equation (12) will become 

1 gsin(^-Q) ) 


1 

■a(l-e': 


0 1 l/l + 7* iin*(v"^Ti') 


and this equation gives 

maximum value of r = a(l — e®) 




minimum value of r = a(l — e® -i- J 1 + 


It will be easily seen that the maxima and minima of r derived firom equa- 
tion (13) differ from the true elliptical maxima and minima by quantities of the 
third order. 

Let us now put the first differential coefficient of equation (12) equal to 
nothing, and we shall find, 

sin (u — 0 ) ^ / cos {v — (o) cos (v — Q) sin (v—Sl) _ q 

7^dv y^l +^8m®(v'^3) {l-h7®8in®(v — g3)}i 

It is evident that this equation cannot be satisfied when v — <u = 0 unless 
v — Sl also equal 0° or 90°, which requires that- the transverse axis of the orbit 
should be in the line of the nodes or in a line perpendicular to it; and as this is 
not the case in nature, it follows that the equation is not applicable to the problem. 
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10 Let us now transform the equation of the ellipse in such manner that the 
anomalies, instead of being measured on the plxne of the orbit, may be measured 
on the fixed plane which is inclined to the plane of the orbit by the angle %, and 
let the place of the ascending node of the orbit on the fixed plane be denoted by 
Q, while the quantities v, 6 and cu have the same significance as in the last 
article , also let v' and a' denote the place of the moon and of the perigee m the 
orbit, and 0„ the latitude of the perigee Eeference to the figure will make the 
geometrical conception clear 

The equation of the ellipse will evi- 
dently give 

1 _ 1 + e cos {v' — Q}') 

r a(l-7) ' 



(17) 


and the required transformation consists in finding the value of cos (v' — q>') in 
terms of the angle v — o) Now, if we put v' — a>' = (v' — Q) — (<o' — Q), we shall 
have 

cos (v' — (o') = cos (v'— Q) cos (co' — Si) + sin (v' — Q) sin (o)' — Q) 

We also have 

tan a = y, tan 0 = y sin (v — Q), 

1 


coad = 


cos 05 = 


V'l + y"sin*(v — 53) 
1 


■j/1 -l-y^sin“(i) — £3) 


(18) 

(19) 

( 20 ) 


|/l-by'‘‘sin“(<!»— £3) 


sin 00 = 


y sin ((o — £3) 
•l/H-y*siu(cti— £3) 


( 21 ) 


Then from the right angled spherical triangles of the figure it is easy to 
deduce 

cos (v' — £3) = cos (v — £3) cos 0, sin (v' — £3) = sin 0 cosec ^ , (22) 

cos (a)' — £3) = cos (<o — £3) cos 0o , sin («)' — £3) = sin 0o cosec z (23) 

If we substitute these values in equation (18), it will become 
cos (v' — (o') = cos (v — £3) cos (a> — £3) cos 0 cos 0o + sin 0 sin 0o cosec® i (24) 


Substituting the values of sin 0o and cos 0o in this equation, it becomes 

, , ,, cos (v — £3) cos ((0 — £3) cos 0 + y sin 0 sm ((o — £3) cosec® i 

cos [v' — (o') = — ^ 

l/l + y® sm® (a> — £3) 

Now we have y sin 6 =y^ cos 0 sin (v — £3), 

and y®, cosec® ^ = tan® i cosec® i = 1 + tan® « = 1 + y® 


(25) 

(26) 
(27) 
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If we now substitute the value of ysin(? in equation (25), and then the 
value of "f cosec® i, we shall obtain, after multiplying by e and putting 

cos !?(, , 

l/l + -f sin® (<0 — Q) 

e cos (a' — w') = e cos (?o cos cos (u — gj) cos (« — Q) 

+ e (1 + -f) cos 00 cos 0 sin (u — Si) sin (<t» — S^), . . ^28) 

= e (1 + If) cos 00 COS0 cos (v — w) — |ey® cos 0o cos 0 cos (?; +_w — 2 £3). . 

Substituting this value of e cos {v' — m') in equation (17), it becomes 
11 

- = -j-jj { 1 + e(H- ^f) cos 00 cos 0 003 ( 1 ;— fu) — •le)'® cos0o cos0 cos('e+«; — 2,£J) } , (29) 

which is the polar equation of the ellipse, in which the polar angle is measured 
on a plane inclined by an angle whoso tangent is y, to the plane of the ellipse 
itself. 


Lot us now discuss this equation and see if it meets all the requirements of the 
problem. If wo substitute the values of cos0 and cos0o, in the second member 
it will become 


1 

r 


a(i^7){^ + 


__ e cos (v — m) — I ef cos (v -f <0 — 2,Q) 

1 / 1 -1- f sin * (aji —~Si) \/l +f sih^^iii -- S) 


(30) 


If we suppose that v ~ to, it will givo 

.i-^J_«>'*” i«>'*co.s®(w-£5) , <l-I-y®sin®(ft;-a) , , 

whatever be the relative values of to and Si. 

The maximum and minimum values of r are also equal to a(l + c) and ffl(l - e), 
which are their correct values. If wc now take the differential coefficient of 
equation (30), wo shall find that it may be reduced to the following form : 

fooN 

r^dv. {l4 y*Hin®(i;--Q)} i ’ 

This always becomes equal to nothing when v = to, and when v ~~ to + 180°, as 
it manifestly should, since the radius vector is a maximum or minimum at the 
extremities of the transverse axis. 

If we now divide equation (29) by cos 0 « — --..i it will become 

■j/l + y®8in* (u — £'J) 
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1 I {-/l + f sin'C'i? — + if) ■“ 

rcosd a(l— e*) 

— ^ey^cos^ocos (v + aj — 2S2)}, 

This IS the correct value of u corresponding to an elliptical orbit, and if we 
suppose that ai = ^3 (m which case cos dt = 1), it will become 




u = — = — {T/l + y^sin'‘(i;— &3) + ecos (v — o)}, (34) 

a(l— e“) 

which IS the same as eq.uation (3), after substituting the value of i/l + s^ Ec[ua- 

tion (33) gives vTT? 

maximum value of -w = 

a{l — e) 


and 


minimum value of u ■ 


(36) 


a(l + e) ' 

which are identically the same as before found m equations (7) and (8) 

We thus see that equations (30) and (33) meet all the requirements of the 
problem, and are therefore correct For the case in which o) = Q>, they immedi- 
ately change to the equations used by La Place and Plana We also see that 
for the general problem their equation is erroneous by terms of the third order, 
and hence their development of the lunar theory must also be erroneous by terms 
of the third order 


Having thus shown that the equation which determines the moon s distance 
from the earth, according to the theoiies of La Place and Plana, is not a conic 
section, except for a particular case, which does not exist in nature, we shall now 
deduce from the general differential equations of the motion of a body acted on 
by the forces of gravitation, the true law of its motion when the positions and 
magnitudes of the different forces are given 



CHAPTEE I. 


GENERAL DIFFERENTIAL EQUATIONS OF THE MOON’S MOTION, WITH TEE THEORY OF 
HER ELLIPTICAL MOTION, AND THE VARIATION OF THE ARBITRARY CONSTANTS. 

1 . For this purpose let us take the general differential equations of the 
motion of a body which is acted upon by any number of gravitating forces, which 
equations are as follows {Mecaniqm C&este [499], Bowditch’s translation) : 


ddx Id 

IQ\ 

ddy 

/dQ\ 

ddz /( 

cZt® ~\6 

'xj' 

df ~ 

Uy/' 



(A) 


In these equations x, y and z denote the rectangular co-ordinates, and the 
function Q represents all the forces which act upon the body whose motion is 
required. The time is denoted by whose element dt is supposed to be constant. 
We may put equations (A) under a more convenient form for computation in the 
following manner. If we denote the sum of the masses of the moon and earth by 
fx^ and all the other forces which act upon the moon by i2, the function Q will be 
given by the equation 

q = ( 1 ) 

r 

The co-ordinates, x, y and z, of the moon will also he given by the equations 
a: = rcos^ cosv, y = rcos0sinv, z — r&ind. (2) 

If we substitute these values of x, y, z and Q, in equations (A), they will 
become 

ddr — rdv^ d — rdff^ , /j. fdB\ ,o\ 

( 8 ) 


2rdrdv cos* 6 — 2r* sin d cos d dvdd + r* cos® d ddv 


di? 

2rdrdd + r* sin d cos d dv^ + i^ddd 

df 


Uuj' 


UW’ 


( 4 ) 

( 6 ) 


(AO 


The integrals of these equations will be the polar co-ordinates r, v and Q of 
the body whose motion is required. 

- If we multiply equation (4) by dt, and integrate, we shall find 


r® cos® 




( 6 ) 


c being an arbitrary constatit quantity. 
3 
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If we now multiply equation (4) by tan^ sin-udi, and equation (5) by 
cos^d^, and take the sum of their products, we shall get, after putting 
2 sin^ 0=^1— cos^ d + sin^ 6, 


^ f 2rdrdv sin d cos ^ sin v + r^dv^ sin 6 cos 5 cos v + r^ddv sin 6 cos 0 sin-u ' 
^ ^ + 'if^dvd d cos^ ^ sin v + 2rdTdd cos v — r^dvdd sin v + T^ddd cos v 
I — r^dvdd s^m^d s,mv 


-c?4(— 


tan 6 am V + 


(dB\ ] 

\d&) j 


(7) 


Equation (7) gives by integration. 


— {r^dv sin 0 cos ^ sin y + cos y } 

dt 




( 8 ) 


In like manner, if we multiply equation (4) by — tan d cos vdt, and (6) by 
sin vdtj we shall get, by integrating the sum of their products, 


dt 


\fd6 siny — r^dv sin d cos 0 cos y } 


( 9 ) 


e' and c" being arbitrary constant quantities 

If we now multiply equations (A') by r^{cZy cos^siny + cZ^sin^cosy}, 

_|_ 2rcZy cos ^ cosy, and sin ^ cosy + 2cZ^ cos ^ cosy, respectively, and 

cos^ 

take the sum of the products, we shall obtain 


' Sr^drdd^ cos ^ cos y + 2'^d0ddd cos 0 cos v — r^d0^ sin 0 cosy 
— 7^d0^dv cos 0 sm y + Sr^drdv^ cos^ 0 cos y + 2r^dvddv coa^0 cos y 

1 —3T^dv^d0 sin 0 cos® 5 cos y — r^dv^ cos® ^ sin y + 2r(ir®dy coa0 sin v 

J y 

fZi* 1 + r^ddrdv cos 0 sm ® + i^drddv cos ^ sin u — i^drdvdO sm ^ sin v 
+ T^drdv^ cos ^ cos u + ^rd-r^dd sm ^ 0081 ? + i^ddrd d sin 6 cos v 
+r^drdddsmdcoav + 7^drdd^ cosdGOSv — 'i^drdvddamdBinv 
+ fi{coa 0 sin vdv + sin 6 cos vdd) 


— —i^{dv Goad smu+cZ^ sin^ cosd} +2rdv GoaO cosu} 

\dr J cos 6 

{dr sm ^ cos t; + ^dO cos 0 cos v} 



( 10 ) 


Multiplying equations (A') respectively by r*{£Z^ sin <? sm u — cZi) cos 0 cost)}. 
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— + 2Tdv cos 6 sin v, and dr sin ^ sin + 2Tdd cos 6 sin v, the sum of their 

cos 0 

products will give 

Zr^drcW^ cos ^ sin v + 2T^d6ddd cos /? sin v — r^dO^ sin d sin v ' 

+ rhhdff^ cos d cosz? -f ZiMrdv^cos^O sin^y + 2'tMvddv cos^d sin?; 

1 — ZrHvHO si nd cos^d sin 4- r^dv^ cos^O cos v — 2Tdr^dv cos d cos v 

^^2 ' — rhldrdv cos 6 cos v + rdrdvdO sin 6 cos v + cos 0 sin v 

— 'if^drddv cos 0 cos v + 2Td/)hW sin sin ?; + rhldrd 0 sin 0 sin v 
+ rhlrddd sin 0 sin v + T^drdO^ cos sin ?; + T^drdvd0 sin cos ?; . . (11) 

+ fjL{d0 sin 0 sin — dv cos 0 cost’} 


■T^ldvcosdcosv — dO sin/? Sint?} 


+ { dr—— — 2rdv cos 0 sin v} 

cos 0 \dv 


‘ {dr sin 0 sin v + 2TdO cos /? sin?; } 


If we now multiply equations (A') respectively by — ?’^6Z/?cosd, 2TdvsinO, 
and — dr cos 0 + 2Td0 sin /?, the sum of their products will give 

C — cos 0 — r^dd.nl0 cos 0 — i^dixldO cos 0 -h sin /? "^ 

^ J + 3A/';Y?/9-sin 0 + rhhhW cos*’/? + 2'ii^dMdd0 sin ^ + tH0^ cos 0 I 
^ 4' ZrHnlv^ sinO cos^0 4- 2T^dvdd/D sin 0 cos^O — 2iHvH0 sin®0 cos 0 [ 

-fxco^odo r 


= 7 ^df)CQsO 


— 2 rdvhm 0 


+ {dr cos B — 2rd0 siu B} 


Lastly, if we multiply ct|iiatiou8 (A') respectively by 2dr, 2dv and 2dd, the 
sum of their products will give 


2drddr + 2rdrdv^ cos^d -I- 2rdrd(P + 2'r^dOddO 
+ 2't^dvddv cos* 0 ~ 2i^dv\Ui siu 0 cos 0 


+ 2/4 


Kf)^- 


-(fh 


If we now take tho intogrnls of equations (10), (11), (12) and (13), we shall obtain 
-j{r*c^^*cosdcoBv+ry'!;*eo8®dcoa'y +r®fW'»coB5Bini? +r*cZj*cZ^ Bind cos v} 1 
— /tcos^ cosv—/ J 

■ — J' |r*{cl'ycosdsin'!;+cZd8iiidcosu}(^^}+(cZr“™+2rcZucos0cos^)} 

+ {dr Bind COS'D + 2rdd cosd cosu} | 
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i - {r®tZ5*cos^sm'u + r®d'y®cos®0smv — r^drdv cos 0 gosv +r®cZrcZ(?sin0smi;} ^ 
d^ 

— fi cos 6 smv —f' 

►— ^ \'t^{dd sin0 smv—dv cosd cosu} |'^^+{2rc£t; cos^ sirw)— 


J I ' 


'\dr} 


+ {dr smd Binv + '^dd cos 0 smv} | 


( 16 ) 


de 


{r^dd^ sm d + r^dv^ sm d cos* 6 — r^drdd cos 0} — jU sm ^ —f' 


=^|r*(i0cos^ ~ 2i*d!vsin0 {^^^—{^dd smd — dr cosO]i^^ | 


, ( 16 ) 


+ rW + r*(^v* cos* fl} - ^ + 


a 


= -2 




(IV) 


f) /^ arbitrary constant quantities to complete the integrals 

Equations (6), (8), (9), (14), (16), (16) and (17) are the polar equivalents of 
equations (P), Meeamque (Meste [672] 

If we now multiply equation (8) by cos v, and equation (9) by sm v, the sum 
of their products wiU give 


dd e'cosv + o^^smv 
dt r* 

Equation (6) gives 

^ _ __c 

dt r* cos* 0 r* cos* dJ\dv) 


(18) 


(19) 


If we now multiply equations (14) and (15) by sinv and — cosv respectively, 

( 3/0 

and divide the sum of their products by r* cos d — , we shall get 

dt 
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~ _ /sini) —f'cQS.v dt 
dt r^cosd ‘ dv 


r®cos0 

r 

COS vdt . 

r. 




C< 


T^Goaddv 


7^{dd sin ^ sin V — dv cos 6 cosu} 

+ 1 2rdv cos 5 sin® — dr^^ ]■ \ 

1 cose i\dvj 

+ {2rdd cos 0 sin u + dr sin 0 sin ®} 

7^{dd sin ^ cos® + dv cos^ sin®}^^^ 

1 2rcZ® cos 0 cos ® + ]■ / 

1 cose J \cZ® / 

+ {2rcZ^?cos^cos® + cZrsin^ C03®}^~— ^ 


( 20 ) 


The integrals of equations (18), (19) and (20) will give the latitude, longitude 
and radius vector of the moon at any time t. To determine the integrals of these 
equations it is however necessary to know the values of the constant quantities 
e, <t', a", f, /', /" and a which were introduced by the integrations, when the 
function JR is equal to nothing, and we shall now proceed to determine them. 

2. For this purpose we shall suppose that JR 0, which gives — 0, 
0, and 0; and equations (6), (8), (9), (14), (15), (16) and (17) 


will become 


r®cos®0™ 1 = c; 
dt 

r® sin 6 cos 6 sin ®” + r® oob®~ : 

dt dt 


■o'; 


dd 


dv 


r® sin®~ ~ r* sin ^ cos 9 cos® . 
dt dt 


_1 f cos ^ cos ® + r*d!®* cos® I? cos® I „ _ 

ciIt^l + r®tfr£Z®cos(?8in® + r®drcZ^sin^co8®J •’ 

_1 Jy^cZd^^cosdsin® + r®cZ®*co8*^Bin» 1 „ • —ft. 

1 — r^drdv cos ^ cos ® + t^drdd sin ^ sin ® J ^ ^ ^ ’ 

{r’cfd® sin d + r®(i®® sin 0 cos® 6 — r^drdd cos 9} — fi sin 9 =f " ; 


( 21 ) 

( 22 ) 

(23) 

(24) 

(26) 

(26) 
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{cZr* + rW + cos* ^} - ^ + ^ = 0 


Equations (18), (19) and (20) xlso become 

dd g'cosv + e'^smu 
dt^ r* 

dv 0 


dt r*cos*^ 

dr fsmv —f' cos v dt 
dt 7^ cos d dv 


(27) 

(28) 

(29) 

(30) 


If we substitute the value of dt given by (29) in equation (28) it will become 


dd 

cos*0 


= — cos vdv + — Quivdv 


Equation (31) gives by integration 


tan 0 == — sm V cos v 

c c 


(31) 


(32) 


When the moon is at the node v — Si, and 6 = 0, therefore equation (32) 
will give 

At ! 

(33) 


e' e" 

0 = — sinSS — — cosQ, 
c c 


tlijsrefore 


0 

tan Q> = — 


(34) 


dd 


When ^ IS a maximum it is equal to the inclination of the orbit , and -j- = 0 , 


equations (22) and (23) will therefore give, by substituting the value 


dt’ 


dt 


c' = c tan d sinu, d' = — c tan 6 coav (35) 

If we square these equations the sum of their squares will give 


tan*^ = 


c'* + e"* 
c* 


= f 


(36) 


The quantities e, c' and c" therefore determine the place of the node and the 
inclination of the orbit , and if we now put 


■= ycosfi^. 


— = ysmQ, 
e 


(37) 
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eq[iiatioii (32) will become 

tan 6 = y {sin d cos — cos v sin Q} = y sin {v — Q), 

•which IS the same as equation (2), Int , and gives 

^ y sin (d — £^) o 1 


sinp = 


l/l + )^sui^(v 


l/l+/sin’‘('U — S2) 


(38) 

(39) 


If "we now square the values of e, c' and e", and put their sum equal to A®, "we 
shall find 


This equation gives 


j/2 ^ Qin — cos^ 6<M + 

at/ 

(40) 

cos* h? 

df ~ r** 

(41) 

equation (27) it •will become 


df 7^ r a 

(42) 


At the extremities of the transverse axis, dr = 0, and equation (42) wiU 


give 


This equation gives 


and 


fxr^ — 2/iar = 

= ah 

(43) 

T =a + a\i 

1^. 

(44) 


jua 


r = a~a\j 


(45) 


f/a 



The sum of these values of r is equal to the transverse axis of the orbit, and 
their difference is the double of the eccentricity Therefore 2a denotes the trans- 
verse axis of the oibit, and -k/l — — denotes the ratio of the eccentiicity to the 


semitransverse axis or mean distance If we put 


' /ua 

(46) 

equations (44) and (46) •will give 


r — a(l ±:e), 

(47) 


which are the values of t at the extremities of the transverse axis 
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If we now denote the longitude of the moon, when is a minimum, by <y, and 
also put d/f 0, and v — <Oj in ecjuations (24) and (25), they will become 

■— cos d + r W cos®d} — fi cos d j cos o>=f, (48) 

~ cos d + cos®d} — ju cos d\ sin a)=f (49) 


These two equations give 


tanai = 


/' 

/ 


(60) 


If we now take the sum of the squares of the values of /, /' and f" given 
by equations (24—26), we shall obtain 


+ d't^ coa^ey + rW(dd* + dv^ cos"(?)} ] 
air 

But equation (41) gives 


( 61 ) 


d8^ + dv^cos‘d = 

Substituting tbis in equation. (61) it becomes 


(52) 


(63) 


If in this equation we put dr = 0, r = a(l± e), and h? = au(l- e*), it wiU 
become 

/*+/'* = (54) 

If we now substitute r = a{l-e), v = (o, 6^8, and«Zi* = 0, 

0 

in equations (24), (25) and (26), they will become 


r 

dO^ 

1 

1— 1 

[ COS^^ocZu* 

f ^ 

dd^ 

Lo(l — e) 

[cos^tfod'U* 

«■ ! 

f dd^ 

Lo(l-e)l 

1. cos* 6^v* 


1 ■ 
COS*^o, 

//jcosdocosa* 


(55) 

1 1 
cos“^*J 

j- — yujcosdosineo 


(56) 

-i-1 

cos® do J 

■ — ;u|smdo =/" 


(57) 
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Bub equation (38) gives, when = 



dd = y cos® do cos {o> — Q) dv, 

(58) 

m 

whence we get 


(69) 

and 


(60) 

whence we get 

+1 ! + »■ 

coa*dodv^ cos® do 

(61) 

and equations (66-67) become 



f ^(1 -h _ 1 cos cos fi> =/, 
(.a(l-e) J 

(62) 


1 ^ 1 cos do sin o) = /', 

(a(l-e) J 

(63) 


{4i4 ' 

(64) 


But equations (36), (40) and (46) give 

e“(l+/) = c» + e«+c"® = ;i® = a/4(l-^), (65) 

therefore we get — ~ (®®) 

o(l — e) 


Substituting this ^alue in equations (62-64), we get finally, 

f = fxeooBdoGoso), "I 
/' = fie oos do BID. (i>, y 

f — fxeBmdo. J 


Therefore the quantities /, f' and denote the pDroduct of the sum of the 
masses of the moon and earth into the co-ordinates of the centre of the orbit 
when referred to the focus as the origin. 


3. Having thus found the values of the constant quantities introduced by the 
integrations, if we now substitute them in the differential eqimtions of the co- 
ordinates r, V and d, we shall obtain, by means of another integration, the values 
of these co-ordinates in terms of the time. But as we have already found the 
value of d in terms of v, in equation (38), we shall also find the values of r and t 

4 
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in terms of v "We shall then, by inverting the formulas, be able to find the 
values of r, v and d in terms of the time t 


Eq^uation (66) gives 

l/TTY 


( 68 ) 


If we substitute this value in eq^uation (29) it will become 

dt i/l +-f , , „ 

— i ‘TCOS^d 


(69) 


And if we also substitute these values, and also the values of / and in 
equation (30), it will give 


dr e(l + y^coso to , \j 

— = ■ ■ „ cos®^ sin (v — 0)) dv 

r a(l — e‘) 


But we have 


cos® a = 


{1 + y®sin®(i; — ^i)}^ ’ 


(70) 

(71) 


which, being substituted in equation (70), gives 

^ _ e (1 + y®) cos ^0 sin (v — w) dv 
r® a(l — e®){lH-y®sin®('y — SS)}^ 


(72) 


This IS identically the same as equation (39), Int , and hence equation (72) 
will give, by integration. 


1 

r 


a(l 




1 + 


e (1 + ^y®) cos dp cos (v — w)— cos d^, cos (i? + 6) — 2$^) 
■j/l + y* sin'* {v — SS) 



which IS the same as equation (36), Int , obtained from purely geometrical con- 
siderations, and IS therefore the equation of an ellipse 

If we now substitute this value of r in equation (69) it will become 


^ a^(l-e®)Vl-t-/co3®g 

i/Ti[ 1 -1- ^ ^ ~ ~ ^ “h ~2^) I * 

^ I i/l + y®sin'*(u — f2) J 


(74) 


This equation expresses the rigorous relation which exists between the element 
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of time dt and tlie differential of the moon s motion m longitude corresponding 
to any part of the orbit , and its integral will give the true time t, which is 
required for the moon to pass through any arc of longitude which is denoted by v 
It IS not, however, readily integrated in its present form, but we may develop the 
second member into an infinite series arranged according to the ascending powers 
of e and y, which, in the lunar theory, are numerically small quantities, and carry 
the approximations to any degree of accuracy which may be necessary The 
different terms of the equation when thus developed can then be integrated 
separately, without any analytical difficulty We shall now attend to this trans 
formation of equation ( 74 ), and shall carry the approximation to terms of the 
seventh order of magnitude depending on the eccentricity and inclination of 
the orbit This degree of approximation is greater than is necessary in the 
theory of the moon a motion , but as the same formula may be applied to the 
motions of the planets, it was thought best to give it all needful extension for 
that purpose 

4 If we develop the variable part of the second member of equation ( 74 ) by 
the binomial theorem, it will become as follows 

q 4. 6 (1 + jry^) cos 60 cos (u — w) — cos dp cos (v + o — 2Q) 

•j/l + y* sin'^ (u — iJ) 

os* ^ — 2 ('^ ~ ~ (i) + a> — 2Q,) 

{1 + y*sin*(u — gS)}i 

Q {e (1 + j^y*) cos ffo cos (v — w) — ^ey*cos ff p cos (v-hw — 2Q)}* 

{1 +y*siu*('y — S^)}* 

^ {e (1 + ^y*) cos dp cos (v—w)— ^ey*co3 dp cos (1; + a) — 2Q)}* 

{1 + y*sin (y — 53)}^ 

K {e (1 + ^y*) cos do cos (v — to)— ^ey* cos Op cos (v + w — 2 g 3 ) }* 

{l+y*8in (u — 53)}® 

_ 0 {g (1 + ly*) cos dp cos (v — w) — ^ey* cos ffo cos (v + — 253 ) 

{1 +y*sin®('y — 53 )}^ 

17 {e(l + ^y*)cos(?ocos('y — tt>) cos (i> + a> — 253 )}® 

{l+y*sm*('u — 53 )}* 

— 8 {g (^ "I" go** ^0 cos (1’ ~ <»>) ~ coso cos (y + — 253 ) }^ 

{H-y*sin*(y- 53 )}^ 
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We shall now develop each term of this equation separately Since 

ooad = ^ we shall obtain 

1 / 1 + 7* sin* (t) — £3) 

cos* ^={l+y*sin*(« — £S)}“‘ = 1— 7* sin* (v — £1) + y* sin* (v — £1) — / sin* (v — £^ ) 

= l-i7*(l-|y* + |7*)+i7*(l-y*+if/)cos2(i;-£3) ') 

+ i/(l~|/)cos4(«) — £5) +-^y*cos6('u — £S) J ' 


We also have 

C08ffo= {1 +y*sin*(a> — £3)}-i 

= i -if + A/- Mf + if a - if + i^^f) cos 2 (n> - Q) 

+ A/(l“|-}^cos4(tu — £2)+7^7 ®cos6(<u — £3) / 

We shall therefore obtain the following values of the different terms of the 
second member of equation (75) 


^ ~b if) cos &o cos (v — w)— cos^q cos (^? + a> — 2£S) 

{1 + 7 *sin*(v — £3)}J 


— ^0-—if+ liy* — ^f) cos (v—<i>)+ -^ey* (1 — |y*) cos Z{v — a>) 

+ fey*(l-y*+||y*) cos(3v -<o-2Q)+^ (1- coa(i;-3a>+2£3) 

~ (1 ~ -fy^) cos (3u H- tti — 4£3) + ~ if) cos (u + 3<w — 4£3) 

+ A«y^ (1 - if) cos {5v-(o~ 4£3) + (1 - f y*) cos (t> - 6a> + 4£3) 

-7fjey*cos(5v+<u-6£3)+^6y*cos(v+6fl»-6£S)+Tf|,^cos(7u-«-6£3) 
+ TM^ey* cos (5v -3<o- 2Q,) + ^^ey'cos (3-0 -5a> + 2Si) 

+ TTsVr^y' cos (v—7(o + 60,) 


, ( 78 ) 


{e(l +^y*) cos (v - m) -|ey*cos (v + n» -2£^)}*cos*^o{l + y*sm* (v - £3)}“* 

= Je* (1 — •ly* + |y*) + |e* (1 — Jy*+ Jy^ cos 2 (v—co) +^e*y*(l— y*) 00^2(1? — O ) ' 
+ ie*y* (1 - y*) cos ( 2 v - 4 a» + 20 ) + -^6*7* cos ( 2 i; - 6n> + 4 £ 2 ) 

+ iVe*y*cos 4 (d - o)) + cos (61; - 2 tti - 40 ) 

+iey(l—/)cos(4v — 2(m-20) 
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{e(l +|-y®) coa{v — to) — Je>'’*cos (v + <o — 2Q)}®cos“^o{l +)'*sin®(t; — Q)}“^' 
= 1^ (1 — + }-y*) cos (v — to) + (1 ~ Jy® + 3 (u — a») 

+ ■jViT^'V °°® 5('W — <i>) + -^eV (1 — y®) cos (2v — to — 2Si) 

-^e'Y cos (3^; + w - 4S^) + (1 - f) cos (^; - 3w + 2g3) 

+ (1 ~ f) oos (5v — 2 to — 2Q) + *(1 — f) cos (Su — 6a> + 2S2) . 

— T6Tf cc® (Su + w — 4S2) + cos (u + 3 w — 4^S) 

+ 3 -ih-®'V'‘oos {v — 5<o + 4S2) + cos {2v—7fo +4^S) 

+ ^®i''je*y*cos {7nt — 2to — 4£3) 


;(80) 


{e (1 + ■Jy') cos {v — co)— ■^cy^cos (t; + a) — 2g2)}‘cos*0o{l + y^sin (v — £J)}“® 
== le* (1 — Jy®) ^- •^e'* (1 — -Jy®) cos 2(v ~ to) + ■J-e* (1 — -ly®) cos 4 (■» — w) 

+ -i^g-ey cos 2(v ~ Q)+ • 3 %sy cos(2u — 4ft» + 2Q,) -l-f^ey cos(4v — 2<u — 2Q) 
+ •S^ey cos (6i; — 4a> — 2^';^) + fg-ey cos (4v — 6ai + 2Q} 


;( 81 ) 


{c (1 + J-y®) cos ('U —to) — |ey®cos (® + w — 2S^)}"cos <?o®{l + y® sin® (^; — Q) }~^ 

= |e® (1 — |y®) cos (v — to) +■ (1 — ^-y®) cos 3 (v — to) 

+ (1 ” J-y®) cos 5 (v — to) + -f^j^g-e^y® cos (3'i; — to — 2£i) 

-1- xIt <^y* cos (■« ~ 3 w + 2S^) + j-l^ey cos (5v —3co — 2Q) ■ - 
+ ■a®2c'V®cos (3v — 5 to +2.Q) + -ji-ffty cos (7nt—Bio — 2Q,) 

+ cos (5u — 7<u + 2£I3) 


(82) 


{e(l H- ^y®) cos (v — to ) — ly wa {v + w --2Q)}*cos”(9o 
{H-y®sin®('W — 

— ■j^y* 4 Jf c^cos 2(v — io) + cos i{v — to) + -^c'cos 6 (■!; — «>) 


(83) 


{e(l + ^y ®) cos ('(? -- to) — co.s { v -H <t) — 2 g^) } ® cos® 

{1 4y®sin®('u — Si)}i 

= -|f e® cos (v — to) 4- I J-c® cos 3 {v — to) + cos 5 (v — «>) 4- ■^ 46 ® cos 7 (v — £S) . 


(84) 


If we now multiply equation (78) by —2, (79) by +3, (80) by —4, (81) by 
4- 6, (82) by — 6, (83) by 4- 7, and (84) by — 8 ; and then add the products to 
equation (76), wo shall get the following development of equation (74) : 
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l//i 


^ 

-e®)Vl + )^ 

= 1 -iy* + |e^ + 1/ + A/ + f^e«-fey + 

1 - -ihf + Hey + He® - ffey 1 ^ 

+ W-if + l^ + i/ + He® - |eV*} coa2(v - o) 

- { e* — ^ey + + ^eV — — Hey} cos 3 (u - w) 
+ {fe* — Hey + He^* + He®} cos 4:{v — a>)+ ^e® cos Q(v — co) 

~ {|e® + ie'^ — He®/ + ifeV } cos 5 — <w) — cos 7 {v — a>) 

+ {i-/~i-/+|ey + M/“|ey + Hey} cos2('j; — £S) 

+ i/{l ~ f/} cos 4 (d - Q) + H/oos Q(v — Q) 

- {|e/ - H/ + fey - fey + Hie/ + Hey } cos (3iJ - o - 2g5) 

- {jy — iy + + HW® + M-e®/} cos (v — 3a> + 2gS) 

+ {Hey® — xHe/ + Hey} cos {Bv + w — 4S2) 

“ {Hey® — rHey® + He®y®} cos (v + 3 at — 4Q) 

~ {ifey® — ^e/ + iHey } cos {5v — <o— 4^) 

~ {He/ — xHe/ + tH^/} cos (■!? — Scu + 4^) 

+ Hrey® cos (5v + a)-6Si)— H^y cos (v + 5cu — 6g^) 

- Hiey® cos (7nt - a> - 6Q) - Hrey® cos (u — 7a> + 6fJ) 

+ {|ey — IH® + Hey} cos (2v — 4:<o +2S2) 

+ {|ey — fey + Hey} cos (4v — 2<o — 2Q) 

Hey cos {2v — 6o}+ 4Q) + Hey cos (6v — 2a> — 4S2) 

- {fey - -fey 4- Hrey® + HH*} cos (3 d - 5a> + 2£2) 

- {fey -fey + H^y + Hey} cos (5v-6w- 2Q) 

+ HeV cos (6 d — 4«) - 2g3) + He®/ cos (4nt - 6co + 2gJ) 

~ iHey cos (3n.i5 — 7ai + 4S3) — HVey cos (7v -wBw— 4S3) 

— Hey cos {7nt — 5a» - 2£2) - Hey cos (5v — 7a> + 2Q) 
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Now we have 

+ (86) 
Multiplying equation (85) by this value of (1 — e*)Ji/l + -f, it becomes 
^ ■ £ = 1 “ 2e { 1 — + tI ^ey + ^y«} cos (« - aj) 

+ — ^y*} cos 2{v — <o) 

— {e® + |e* + ^ey* + — -^^ 6 * 7 * — -^ey*} cos 3 (u — a>) 

+ {fe* + + fe®} cos 4 (u — <u) + -^^cos 6 (u — u») 

—{16®+ + -J^eV} cos 5(v — a)) — ■Je'^cos 7 (v — a>) 

+ 1 “ ly® ~ fe* + ^y* — h* + V} OOS 2 (d — £5) 

+ -J-y^{l — / — fe®} cos4(y — Q) + ^y* cos 6 (u — £^) 

— f ey*{ 1 — — -J-e® + + ^eY) cos (3u — a» — 2Q) 

— Jey® {1 - ^y® - ^e® + cos (i; - + 2Q) 

»■ • (° * ) 

+ A®7^{l~7*“■i■®*}co3(3^J + a>— 4£5)— gij-ey^{l— y®— ^e®}cos('y+3a»— 4£3) 

— w — 4Q) + ■ 5 ^ey® cos (5v + o — 6 £^) 

+ 4£3) — ■j^ey® cos (r + 6 cy — 6 £ 5 ) 

+ |6®y® { 1 — iy® — -Jre®} cos (2v — 4<d + 2£3) — cos (7v — w — 6£3) 

+ |e®y® {1 — ly® — Je®} cos (4u — 2(o — 2Q) — •s-f^cy® cos (■« — 7w + 6g3) 

— {|e*y® — ■ 3 ^e®y* + ^-fij-ey®} cos (3v —5(i> + 2Q) + ■ 3 ^ 6 ®y® cos (2v — 6 a> + 4£3) 

— {|ey — -j^ey + ^i'^jj-ey®} cos (5v — Sco — 2Q) + ^eV cos (6v — 2af — 4Q) 

~ cos (3-!; — 7o + 4£3) — ■^■^eycos (7v — 3<w — 4£J) 

+ cos ( 6 u — 4<o — 2Q) + -j^ey cos (4® — 6 (t» + 2 Q) 
cos (7v — 6 <y — 2£3) — .J^ey cos (5v 7<y + 2 Q) 

If we now put n = (88) 

* ■ 

equation (87) will give, by integration, 
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nt = v -2e(l + xf^ey + sm (y — co) 


+ fe® (1+ -K + - i/) sm 2 (v — a)) 

— (K + ie ® + - -h^) Bin 3 (u - a)) 

+ {^«* + -h^ + A«V} Bin 4 (« — a)) + y^e'sin 6 (i; - a>) 

“ { Ab® + ^ + - 5 ^*} sin 5 (v — a>) — sin 7 (u — a>) 

+ iy^U “ ~ + l®V} sin 2 (u - ft) 

-i/{l-i/“l^ + T^/“i«‘ + t®V}ain2 (a>-g2) 

+ { 1 - / - 10*} sm 4 (u - SS) + xk/ sm &{v-Qi) 

~ •sV/ { 1 - / — 1«®} sm 4 (a» — Q) — xix/ sm 6 (a> — g2) 

-K{1 + w - + isV} sm (3v - a, - 2gJ) 

— i6}^{l — i-/ — -K + H-/ “ As* + is®7*} sm (u — 3a) +2gJ) 

+ ■g'^ey*!!— y*— ^e*}sm(3D + a» — 4i2) — y*— J^}sm(v+3a>— 4£2) 
~ y® ■" «*} sm (S'!) — a) — 4£J) + sm (5i) + a) — 6SS) 

— Asy®{l — y* ~ e*} sm (v — 5(0+ 4gJ) — (v + So) — 6£S) 

+ '^6V{1 — -Jy® — Je“}sm(2'!)— 4a)H-2S5) — xfxsy'sm (7 d — a) — 6SS) 

+ AsV{ 1 ~ iy* ~ i^} sm (4 d — 2a) — 2gS) — -B-fysy® sm (v — 7a) + 6gJ) 

~ {i^V ~ TA®y®}sm (S'!) —5a) + 2S3) + AA® sm (2'!) — 6a) +4S2) 

~{is®y*“AsV^ + 'BAsy®}sm(5'!)- 3a)— 2S3) + A«ysm (6 d— 2a)— 4£3) 
— xf?e®y* sm (3i) — 7a) + 4Q) — xA^® sm {7v — Sco — 4£S) 

+ A®V sm (6'!) —4a) — 2g3) + A®V sm (4-1) — 6a) + 2Q) 

— A^V sm (7'!) —6a) — 2Q) — A®V* sm (6i; — 7a) + 2Q) 


(89) 


In this equation nt denotes the mean longitude of the moon, and v the true 
longitude , and the con'itant introduced by the integrations has been made to 
satisfy the condition that the mean and true longitudes shall be equal to each 
other at the extremities of the transverse axis By inverting this series ’we may 
obtain the value of v m terms of nt 

6 For this purpose we shall observe that if we change the sign of all the 
terms of the second member except the first, and at the same time change v into 



DEVELOPMENT OE THE ELLIPTICAL MOTION. 17 

nt in those terms, and calling the sum of the terms thus changed /(nit), we shall 
have, according to the theorem of La Ghangb, 


V 


JL I I 1 dflnty , T dj{ntY , i dy(ni 
.. nt +f(nt) + ^ 

, dyjnty d^f jnty d<^f (nty 


(90) 


We shall therefore have 


f(nt) = 2e { 1 - -^y* + + -8^*} («* “ «>) 

— |c®{l + le® + -^e* — sin 2 (nt — (o) 

+ {ie* + ie* + i^ey* + - -^cY — -^ey®} sin 3 (nt — o>) 

~ + tV®* "h ^ ~ ~ 6 (ni — (o) 

+ sin 5 (nt — w)+ sin^ (ni — co) 

— ■Jy® { 1 — iy® — |e* + ~ i®* + 1®^*} 2 (ni! — g2) 

+ iy* { 1 “ iy® — |e* + — \6*‘ H- fey } sin 2 (w — P,) 

— iY { 1 ” y* “ f e*} sin 4 (ni — g5) — ^y® sin 6 (nt — gS) 

+ Ay* { 1 — y® — -fe®} sin 4 (<u — S^) + xA/ si^^ 6 (" ~ 

+ Jey® {1 — ^y® — •J'S® + |-Jy® — + ^Y} si*^ — <o — 2£i) 

+ Jey® { 1 - Jy* - ^e* + fjy^ - ^e* + W) sin (nt ~ So + 20.) 

— -s^ny* { 1 ™ y*— |e®} sin (SjiM- w~4Q) -t- A<*y* { 1— y®— -J e®} sin (nt+ 3tt»— 4^) 
+ As}'* {1 ~ y® — e®} sin (5nt — o> “4&J) — si^' ^ ■“ 

■h y* “ ®®) sin (?ii ~ 5tt> + 4:0) + yAej'* sin («.i -f Gw — 6£2) 

— A*V { 1 ■" 4y® “ ie®} sin (2 tci 5 — 4w + 20) 4- -jey® sin (7nt — 0 — 60) 
” A**y® { 1 — iy® ~ i®*} ?in (4nt — 2w — 20) + -cf-jey® sin (nt — 7w + 6gS) 
+ {-JA®- AA* "i^Tfr^y*} sin(3nJ— 6w + 2Q) — ■Ae®y*8in(2»f— 6w+4S3) 
+ { + 2S3) — AA'*sin(6?n!— 2a»— 4^5) 

4 sin (Snt — 7«) 4 4Q) 4 rl^ffy sin (7ni — 3 w — 40) 

— A«*y® sin (Gnt — 40 — 20) — AeV sin (4nt — 6w 4 20) 

4 AA® sin (7»ii — 6w — 2S3) 4 AA® sin (5nt — 7w 4 2Q) 


(91) 
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Prom tills equation we get tie following values of the different powers of 
f{nt), which are all correct to terms of the seventh order 

f{nty = 

— {|^ + ^ + + •jV’' - — W) cos 

- {2e® — |e" + iV/ - A®V* - A/} cos 2 (wi - a>) 

+ {fe® - liVe® + W - ~ T^\«V “ i/} cos 3 {nt - a) 

- {lie* + + iVey } cos 4 {nt -w)- cos 6 {nt - m) 

+ { A®* "h + AfeA*} cos 5 {ynt — w) + cos 7 {nt — «) 

+ - iieV + tA/} (cos 2 (nt - £S) + cos 2 (o - Q) } 

“ {Ay* — A/ “ AA* } {cos 4 (ni — Si) + cos 4 (<m — £3) } 

— Ary* {cos 6 (nt — fi3) + cos 6 («u — Q) } — xig-ey® cos (39^^ + 3a) — 6S3) 
~‘{e/~iey*+i^ey®“iey “feV + AA*} cos (ni + a) — 2Q) 

+ {iy - i«y* — + AA* — AA’y® + Aey®} cos (nt - 3a» + 2S3) 

+ {ie/ - iy - lY + AA* - -rhY + Aey®} cos (Bnt -a)-2Q) 

r 

+ {iey* “■ iey® ~ iAA*} {cos (nt — 5a) + 4g3) + cos (5nt — to— 4£3) } 
~{ii^“AA*‘“T^*y®+Ti5-y®}{cos(27i«— 4a)+2f3)+cos(4nt— 2a»— 2f2)} 

+ JiiA — liA*“HA*+TAey®5 {cos(3ni— 5a>+2£2)+cos(5TOt— 3a»— 2Q) } 

“ {iey* ” iey' ~ 'i^A*} {cos (Znt + to — 4£3) + cos (nt + 3a) — 4S3) 

+ {Ay* “ Ay® “ A^V} cos (2nt + 2a) — 453) 

— -^izY {cos (6nt — 2a) — 453) + cos (2nt — 6a) + 453) } 

— ifiA* {cos (6nt — 4a) — 253) + cos (Ant — 6a) + 253) } 

+ AVey® {cos (nt — 'J(o + 653) 4- cos (^nt — a) — 653)} 

- xAey® {cos (nt 4- 6w — 653) 4- cos (6nt 4- w — 653)} 

+ in A* {cos (^nt — 3a) — 453) 4- cos (3ni — 7(o+ 453)} 

+ iliA®{cos(7«'i — 5a)— 253) 4-cos(5)it — 7fa4-2S3)} 

+ Apy® {cos 4- 2a) — 653) 4- cos (2«,i 4- 4a) — 653)} 
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/ (nty — {6e® + -l-J-e® + -/jey* + — -fleV } sin (nt — io) 

— {fe‘ + fffe® 4- l-JeV} sin2(ni — w) 

— {2<^ - - U¥—iW “ sin 3(n« - w) 

4 {-fe^ — fje® 4 sin4 (n« - tu) 4 sin 6 { 7 it - cu) 

— {fa-«* - sin 5 (nt -eo)- sin 7 (n« - a>) 

4 { |ey - |eV - -J-l-lcy + tI (t/I {sin 2 (<o — Q) — sin 2 (nt - Q ,) } 

+ {■ 3 %y “ ■sV>'®~ f K/} {sin (Bnt 4 - 4Q) — sin (nt 4 3a> — 4£5)} 

~ ~ f f e*y*} {sin (nt — 5<o 4 4g2) 4 sin (5^^^ — a> — 4^3)} 

+ {^^■e'Y — f J-ffY — ff-gey 4 iVgY} sin (3ni — (o —2Q) 

+ {l-iffy — f 4 -j^/ge/} sin (nt — Z<o 4 2.Q) 

— {f-^c-y — f |ey — e®y* 4 ^-s^ey* } sin (5nt — Zco — 2Q) 

— {fkV - f^y - sin (37 w 5 - Sf«> 4 2g3) 

+ {|e*y» - |ey - + Tf-fiy®} sin (4w« - 2ft) - 2j;3) 

4 {|eV® — |e®y* — + Tfv/} sin (2nt — 4a) 4 2.S3) 

4 ■aVgC®?'* {sin (2nt — Goj 4 4.Q) 4 sin (6n^ ~ 2(0 — 4S3) } 

{sin “ G<o 4 2£3) 4 sin (6nt—4(0 — 2£3)} 

+ { sin 4 (ft) — Q) — sin 4 (nt — Q)} 

“iritry* {sin 6 (« — S3) — sin6 (w! — S3)} 

■i" vl^iry® {sin (2»i 4 4(«) — 6S3) ~ sin (4ni 4 2ft) — 6S3) } 

~ H®VMsin (7n)5 — 3ft) — 4S3) 4 sin (3)ii “ 7« 4 4S3) } 

” { sin (5nt — 7« 4 2S3) 4 sin (7nt — 5ft) — 2S3) } 

+ tSv'V® {sin (5nt 4 ft) — 6S3) — sin (w^ 4 5(o — 6S3) } 

— {sin (7nt — (o — 6S3) 4 sin (ni5 — 7ft» 4 6S3) } 


i 


. ( 93 ) 
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_ {ee® + 3^^ + cos {nt - (o) 

— {86^ — + 1^*} cos 2 (nt — (o) 

+ {96® + IJe^ + H«V} °°® ^ ~ '") 

+ {2e* — + |e®y*} cos 4 (ni — + Jj^e® cos 6 (nt — w) 

_ {3e5 - + HeV} cos 6 (w« - <u) - cos 7 (n« - c>) 

— {6eV® - 3e»)'® + ^eV + ^e/} cos (n,« + a> - 2£S) 

+ {4ey — 2ey - J^eV + T^g-ey®} cos (3«,i — o) - 2£5) 

+ {4eV — 2e®y® — + i^ey®} cos (mi — 3«> + 2£J) 

— {ey — ~ 

— {ey - - 3^ey + -^ey®} cos {&nt-Sa>- 20,) 

+ -le^y® cos (2w< + 2co — 4Q) — -^^y® cos (Snt + 3a> — 6Q) 

— -^ey® {cos (7nt — <» — 60) + cos (ni — 7a> + 6£J)} 

— (cos (5nt— 7ca + 2£S) + cos (7nt— 6a» — 2S3) } 

— ■|^e®y® (cos (Snt + o> — 4S^) + cos (nt + 3<w — 4S^)} 

— {°®® (7nt — Sco — 4S3)} 

+ (2^^ — 6a> + 4SJ) + cos (6nt — 2(o— 4Q)} 

— -^eV {°o® — 4a» + 2SS) + cos (int — 2a> — 2£3)} 

+ (cos (4ni — 6a) + 2SJ) + cos (6nt —Aw— 20)} 

+ -^^eV {cos (nt —6w + AO) + cos (Snt — w — 4£3)} 

+ ^ey® {cos (nt + 5w — 60) + cos (Snt + o — 6Q) } 

— |ey {cos 4 (o) — £3) + cos4 (nt — Q) } 

+ eV {cos 2(w—0) + cos 2 (»S — 
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/ {nty — {20e® + sin (nt — <u) — -^e® sin 2 {nt — m) 

— {lOe® — sin 3 {nt — 0 )) + 16e® cos 4 {nt — <o) 

+ {2e® — + ■|e*y‘‘}sin 5 {nt — <o) — sin 6 {nt — w) 

+ sin 7 {nt — w)+ ^e*y~ {sin 2 (a> — £^) — sin 2 {nt — Q)} 
+ •\S-eV {sin (Sni — 4a> + 2SJ) 4- sin {4cnt — 2w — 2£3) } 

— |eV {sin {int — 6 w 4 2£3) 4 sin {Qnt — 4(t) — 2S2) } 

4 {sin {6nt — 7io 4 2S2) 4 sin (7^^^ — 5(o — 2Q)} 
4X|^(5®J/®{sin(?^^ — 3w 42SS) 4sin(3ni5— <«>— 2SS)} 

— {sin {3nt — 5f« 4 2S^) 4 sin {5nt —3co—2Q)} 

4 {sin (3ni 4 (o — 4S3) — sin {nt 4 3co — 4S5) } 

— {{feV {sin {nt — 5(0 4 4S3) 4 sin {5nt — <o — 4fi2) } 

4 -^g^Y (sin {3nt — 7w 4 4SS) 4 sin {7nt — 3co— 4Q ) } 

/{nty = 20e® — cos {nt — <o) — SOe® cos 2 {nt — to) 

4 cos 3 (jii — (o) 4 126® cos 4 (jii — a>) — cos 5 {nt — to) 

— 2e® cos 6 — 10 ) 4 cos 7 (rii — <u) ~ SOe®}/® C08(rii 4 fa -- 2Sl) 

4 ^6®/ {cos (3ni — <a — 2Q) 4 cos {nt -- 3<a 4 2Q)}' 

— Qe’V® {cos (5n« — 3fu — 2SS) 4 cos (3ni — 5<o 4 2£3)} 

4 l^y {cos (7?1!5 — S<y — 2£3) 4 cos (6n.< — 7<a 4 2Q)} 

/{nty — 70e^ sin {nt — co) — 42c’' sin 3 {nt — «) 

4 146® sin 5 {nt — (o) — 2e^ sin 7 {nt — w) 




22 


THEOKY OF THE MOONS MOTION 


Mnty _ 

ndt 

{ sin {nt - a>) 

+ { 2e® - |e^ + ^Jg-/ - - ^ 4^8 - ^<37®} sin 2 {nt — w) 

- {f «* - • s^«“ + A®/ ~ A®^ “ A®y® “ |-H ®V^} sin 3 (ni - (o) 

+ {fie^ + li®® "I" A® 7 ^} sin 4 (rii — <u) + sin 6 {nt — cS) 

~ {|-2®® + 1 ®^ + sin 5 {nt — 0))— sm 7 {nt — w) 

"" {■H®V ~ tt®V® “■ ■y sV + tA/®} sin 2 (wi — £S) 

+ {Ay* “A/ ~ AA*} sin 4 {nt -Q)+ sin 6 {nt-Si) 

+ {W - i®y* + ■^®y® ~ A®V “ t\®V + A®®y*} sm {nt + w - 2 fi 5 ) 

- {iey® - - ^eV + ^ A* “ tA®®/’’ + TA®y®} sm {nt - 3 w + 2 Sl) 

_ {|®y2 _ |.ey4 — |gy 4. |eSy< _ + -j^ey®} sm {Bnt - to - 2Q) 

+ AV*y®+Tk/} sin(2«< -4^^ + 2Si) 

+ {W- H«V - + A/} sm (47ii - 2co - 2Q) 

~ {flA® — •3^®®y‘ “ A^®®y* + AV®y®} sm {Znt — 5 co ^ 2£S) 

~ { V/®V ~ Hi®®)'* ~ ||®®y® + ■^®/} sm { 5 nt — 3co — 2 Q) 

~ { Ay* ~ A/ ~ A®V*} 8in(27i.i + 2oj — 4S2) + gfgey® sin(3ni + Sen — 6Q) 
~ { A®y* “ A®y® ~ A*} sm (n< — 6a» + 4Q) 

~ { A®y* “ A®y® “ Iff A*} sm {Snt — <o — 4 g 3 ) 

+ { A®y* ~ A®y® ~ •^«'V} sm {Bnt + ft) — 4 Q) 

+ { A®y* “ A®y* ~ tA®V} sm {nt + 3 w — 4 S 3 ) 

+ A? A* sm {6nt — 2(1)— iQ) + -j^A* sm {2nt — 6a> + 4 Q) 

-fHIA* sm { 7 nt -3(0- 4 SS) - A* sm { 3 nt -Y(o + 4 Q,) 

“TllF®y®sin {nt— 7(0 + 6Q) — iAV®y® sm { 7 nt — (o — 6Q) 

FlT®y® sm {nt + 60 — 6fi) + A\®y* sm {Snt + (o — 6Q) 

-iflf® y® sm { 7 nt - 5 (o- 2 Si) - sm { 5 nt - 7 (o + 2 £J) 

“■ Ay® sm {int + 2(0 — 6 S 2 ) — rAy* sm { 2 nt + 4 (o — 6Q) 
ii-i®*y® sm {6nt—4(o — 2 Q) + fff ey sm { 4 ni — 6(o + 2 Q) 
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- {e* + fK + } sin - 0^) 

+ {3e’‘ + -fl-le® + Jf-eV"*} sin 2 {nt — to) 

+ {3e» - We’ + “ A®/ “ ffieV } si^ 3 {nt - a>) 

- {6e‘ — -le® + sin 4 (n.i — w) — W/ e® sin 6 (ni — (o) 

+ { We* - M-M-e' + W/’ey } sin 5 (»< - ^) + ffH^e^ sin 7 (ni -a)) 

+ {W - |ey - fMeV + rhf] sin 2 {nt - O) 

+ f-Jc®)'® sin 4 {nt — Q) — tIt/ sin 6{nt — SS) — AaeV’ sin (2wi — 6fo + 4£5) 

- {H-y- f|y-tAe'y}sin(3nt + <0 -4gS) -^fleV sin(6n«-2<o-4a) 
+ Aey* ~ Aey* ~ AVy} sin {nt + 3^0 — 4Q) 

+ -hey* - Aey® ~ A^seV*} sin {nt - 6fi> + 4^) 

+ {fiy - Mey® - Wi^Y) si“ -a> — iSi) 

- - Aey - AsW + vinY) sin (%< - 3«> 4- 2£5) 

- {f|ey - -ffe-y - fMe®y* + sin (3ni -co- 2 Q,) 

+ { We®y® - We'y “ W/eV + f!iey®} sin {5nt - 3a> - 2Q) 

+ ||( 3 .y _ .g|.eV - + A^ej'®} sin {3nt - 5o> + 2£3) 

- {26®y® — ey — 'W ey “I" A/} sin (4?ii — 2a) — 2£S) 

- {^eV - ieY - m^Y + ihy"} sin {2nt - 4« + 2Q) 

- W ey sin (4ni — 6« + 2fS) — AV ^Y sin {^^ ~ 4a) — 2Q) 

“■ tA/ sin (2n^ + 4a) — 6fJ) + A/ sin (4M.i + 2<o — 6SS) 

+ A® sin (7?ii — 3<() — 4£3) 4- sin (3»i — Yo) + 4S3) 

4- A® sin {5nt — 7(0 + 2Q) 4- eV® sin {7nt — 5w — 2£5) 

- Affey* sin (5nt 4- o) — 6fi3) 4- Arey® sin (ni + 5a) — 6JJ) 

4- H^-ey® sin (7ni — o) — 6£S) 4- sfsY sin {nt — 7(o + 3Q) 
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- {f — ff e® + am 2 (nt — <u) 

+ {-^e® + am 3 (nt — to) 

+ { — -^e® + sm 4 (ni — tu) 

- {-1-1 V — sin 5 (nt — «)) 

+ ^ II ^ e® sin 6 (nt — <o) — sin 7 (nt — to) 

- + -j^jeV + Tl-^ey®} sin (nt+ to — 2 Si) 

+ {f«y ~ f «y — -Vi^sy + sm (Snt — (0 — 2Q,) 

+ {i^y “ — t’^ey H rib®/} s“i — Sw + 2Sl) 

- iW - T^eV - + ^ey®} sin (3nt - 6a> + 2gJ) 

- {Wey-¥W + m^} sm (5nt - 3ca - 2^) 
+ f eV sin (2nt + 2to - 4g3) - -^ey® sm (3nt + 3 w- 6gS) 

~ sin {7nt — a> — 6Q) — yg^j^ey® sin (nt — 7c£» + 6Q) 

- -W^^y sin (5nt - 7«> + 2Q) - (7nt - 6a> - 2Q) 

- fH^V* sm (3wt + <j) - 453) - sm (nt + 3a» - 4Q) 

- Iff eV sm (3nt ~7<o + 453) - sm (7ni - 3<u - 453) 

+ b^y sm (2nt - 6a> + 453) + ^e)>* sm (6nt -2a)~ 453) 

- J^y sm (2nt - 4tt» + 253) - 176V sm (4nt - 2co - 253) 

+ V^y sm (4nt - 6(0 + 253) + if^eV sm (6nt -4^6 - 253) 

+ ■AVey sm (nt - 5a> + 453) + sm (Snt-co- 453) 

- 26y sm 4 (nt - 53) + ffeVsm 2 (nt - 53) 

+ pbey® sm (nt + 5o - 653) + Hf ey® sm (5nt+ to - 653) 
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■ = {16“ + xf + A^Y) sin (nt - <y) - -|e® sin 2 (nt - o)) 

— {-^e® — sin 3 {nt — <t>) + 32e® sin 4 — (o) 

4- + f|-|-eV} sin 5 — w) — ^ e® sin 6 (wi — to) 

+ sin ^ — ft>) — feV sin 2 (nf — S2) 

+ feVsin {2nt — 4w + 2Q) + ^®-eV sin (4ni — 2to — 2Q) 

— |eV* sin (4ni — 6w + 2^) — ^e*y‘ sin (6nt — 4<u — 2Si) 

+ sin (5nt —7(o + 2Q) + gin ( 7 ^^ _ _ 2 Q) 

+ sin (wi — 3<u + 2J5) + sin (3ni — <0 — 2Q) 

— •^eV sin (3»!f — 6<u + 2^) — AJIAe'j/i sin (Bnt —Sa> — 2g3) 

+ sin (3nt + io— 4Q) — xAA* sin (nt + 3(o — 4£5) 

— -g-iTA* sin (nt — 5a> + 4Q) — sin (5ni — w— 4Q) 

+ sin (3ni — 7a> + 4Q) + sin i^nt — 3co — 4^S) 


sin (nt - 0 ) + |e® sin 2(nt — co) 

-’ sin 3 (wi — <o) — •^‘^e® sin 4 (wi — fw) 

+ sin 6 (nt — (o) + •i-fAe® sin 6 (nt — <o) 

+ -^eV sin (ni + w — 2S3) — sin (3ni — co — 2Q) 

— -gVy sin (nt — 3« + 2£3) + -^e®}^ sin (Bnt — Sco — 2Q) 

+ sin (Swi — 6w + 2£i) — ■^|-3^e®y* sin (7nt — 5co — 2S3) 
— -^eV sin (Bnt — 7o> + 2Q,)— sin 7 (nt — a>) 




vrnnr 


n^dt? 


sin (raf — «>) + sin 3 (nt — a>) 

— sin 5 (nt — <o) + sin 7 (nt — a>) 


(103) 


If we now substitute equation (91) and equations (98-103) in equation (90), 
we shall obtain the following value of the true longitude v of the moon in terms 

of the mean longitude nt, 

6 
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V =nt 

+ {2e - -Je® + — ^ + ie/ + sm (nt - m) 

+ {|e® - ^ -ieV -*/} sm 2 (nt - w) 

+ {H e® - i|e“ + W - ¥f + sm 3 (nt - w) 

+ { — tf-K + il«V} sm 4 (n« — at) + sm 6 (n« — at) 

+ { We® - sin 5 (ni-fl.) + sm 7 (?^^ - o>) 

- {if - if -^f + yy +MeV + *7®} sm 2 (nt - SI) 

+ i-hf - A-/ - sm 4 (nt -Si)- ^hf sm 6 (nt - SI) 

+ {if -if- Aey - -sVsV + ^if + AsV} sm 2 (o) - g2) 

+ {■bV/ ■“ -hf — -is^*} sm 4 (<e» — S3 + T^-y® sm &(w — Si) 

+ {ief — ley* — + -i^y* + ief — sm (nt + a> — 2g3) 

- {W - W - + HsV + Aey® - Wyf} sm (3w« -o>-2Si) 

- {Wf - Msy - WsV + */} sm (4n« - 2^u - 2Q) 

- { At/ -f^f-iW) sm (2nt + 2.^ -4Q) 

+ { ~ yyy* ■*■ T%^y } sm (nt — So> + 2S3) 

— {i^y* — iey^ — -IsV} sm (Bnt + to — 4S3) 

+ {iey* — ief — -HeV} sm (nt + 3to — 4^3) 

+ {iy ~ iy ■" W ®V} sm (^nt — to — 4Q) 

— {|-|e®y® — -flsV^ 4- -jig-ey® + W^^V} sm (6nt — 3to — 2£3) 

+ {lAsV sm (2n« — 4to + 2g3) + fJeV sm (Cnt — 2to — 4£3) 

— ey sj n i^nt — 4to — 253) + -^y® sm (4M + 2to — 653) 

_ sm (2nt + 4to — 653) — fi-Hey sm (Znt — 6to + 253) 

— ^ef sm (5nt + to — 653) + Hf sm (nt + 5(0 — 653) 

— -^ey® sm (7nt — to — 653) — sm (Znt + 3to — 653) 

+ sm (7nt— 3to —453) + -^^^ey sm (7nt — 5to — 253) 

+ W ®y sm (5nt —7(0 + 253) + rhe^y* sm (nt — 6to + 453) 
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6. Having found the value of v in terms of nt, if we now substitute it in equa- 
tions (73) and (39), we shall also obtain the values of r and 6 in terms of nt. For 
this purpose we shall first develop equation (73) in a series, carrying the approxi- 
mation to terms of the fifth order, and we shall find 

^ = "i. — e {1 — -^y*} cos{v — to) — i-/} cos (^; •+• — 2S5) 

r 

+ -J-erll — iy“} cos {u — Zto + 2S2) + {1 — cos {Zv — to — 2SS) 

+ tI 8®/ (u — 5a) + 4S3) 4- cos (5u — to — 4S2) 

— xl"sC7'‘ cos (Su -t- o) — 4^) — xf s ey* cos (v + 3a) — 4SS) + cos Z(v— to) 

* 

In order to substitute the value of v in this equation, we shall observe that if 
we put all the terms of the second member of equation (104), except the first, 
equal to jS, we shall have 



V =nt + 

(106) 

This gives 

mv + a = (mnt 4- a) 4- mjS, 

(107) 


m being any positive whole number, and a any angle whatever. 

Equation (107) gives 

cos (mv +a)— cos (mnt + a) cos — sin (7n7i.t + a) sin wij3. (108) 

Since j9 is a small quantity, we may develop sin m)9 and cos m/9 in a series, 
and wo shall have, with sufficient accuracy, 

sin = m)9 — -im”/?, cos = 1 — (109) 

We shall then find, as far as terms of the fourth order, 

/9* = 2e* H- -I- fc’ cos [nt — to) 

— {2e® — cos 2 (nt — to) — 1^’ cos 3 (nt — to) 

— cos 4 (nt — to) — fleV* cos 2 (nt — £5) + ^ 

— .^.jy* cos 4 (nt — Q) — ■^y* cos 4 (to — Si) + -^y* cos (2nt -H 2a) — 4£S) 

— cy^ cos( 7 it -hto — 2£3) I" ^cy^ cos(3nt — o) — ■ 2£3) -1- ■J-ey’ cos(ni—3a)-t-2Q) 

4- — 2a) — 2Q) 4- cos (2nt — 4<<) 4- 2Q) 
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^ = 6e® sin {nt — <d) + e* sm 2 {nt — <u) — 2e® sin 3 (ni — co) 

— • 25 ^e*sia 4 (wf — iu) — fey*sin2(ni — S3) +fey sin2(<y — Q) (HI) 

+ f eV sm (47^^ — 2a» — 2S3) + f«V sm (2^^^ ~ + 2S3) 

^ = 6e^ — 8e‘ cos 2 {jit — o>) + 2e^ cos 4 (ni — to) (H2) 

Therefore eq^uation (108) will become 

cos (mu + a) = cos + a) X 

{ 1 — m® (e* + + -^y*) + — fmV cos {rd, — oi) 

+ (mV — -ImV + ■s’jwV ~ ■frwJ'V) cos 2 (w4 — <u) 

+ l-mV cos 3 (nt — <a>) + (^|^mV H-^jmV) cos 4 (nt — (o) 

+ llm^ey cos 2 (nt — S3) — cos 2 (o — S3) 

+ - 5 *jmy cos4(nt — S3) +-g>jmycos4(a> — S3) 

— < 2 os (2nt + 2a> — 4S3) + \efm^ cos (nt + a* — 2S3) 

— -^m^ey* cos (3nt — <u — 2S3) — cos (nt - 3a> + 2S3) 

— -IJm^ey cos (4nt — 2«> — 2S3) — -^jm^ey cos (2nt — 4a) + 2S3) } 

— sin(mnt + a) X 

{ (2me— Jme®— mV)sin(nt— a)) +(|n!V— ^me*+-jigmy^— f mV)sin 2(nt— a)) 

+ (i-|me® + ^V) sm 3(nt — a)) + + -fniV) sin 4 (nt — aj) 

— (|my® — •|my^ — meV — ■|m®ey) sm 2(nt — S3)— -Jmey® Bm(3nt —w — 2S3) 

+ (^my®— -Imy*— j^meV— ■|m®eV)sm2(a) — S3)+-^ey® sm(nt+a)— 2S3) 

+ ^wy^ sm 4 (nt S3) + sm 4 (a) — S3) — -j^^my* sm (2nt + 2a) — 4S3) 

— (-j^meV + sm (4nt — 2a) — 2S3) — ■Jm®eV sm (2nt — 4a) + 2S3) } 

If m this eq^uation we put in succession m = l, a = — a), m==l, 
a = a>— 2S3, m = l, a= — 3a) + 2S3, and m = 3, a = — a)— 2S3, we aha, 11 
get the followmg eq^uations 
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COS (v — <w) = — e + {1 — -g-e® + ~ ‘it/*) ~ ") 

+ {e — -|e*} cos 2 {jit — «>) + {|-e® — ff-ge'* + cos Z(nt — <w) 

+ |e* cos 4 {nt — w) + cos 5 (»i! — <o) 

- - iW} cos {Znt -w- 25 i) 

+ {iy® “ iy* — -^s^^V} cos (wi + o) — 2S2) 


" { »y“ “ A/ ” ■sW} cos (nt - 3<0 + 2^) 

+ -j-lff/cos (5nt — <0 — 452) — xir/ cos (3?ii{ + <u — 452) 

+ rl^y'* cos (ni + 3 w — 452) — tIt/ cos (nt — 5(0 + 452) 

+ f (jy® cos 2 (nt — 52) — |-ey® cos (4ni — 2(o — 252) 

— ^ey^cos (2nt — 4<i» + 252) — \ef cos 2 (w — 52) 
“ficVcos (5nt — Z(o — 252) — • 5 ® 4 -e®y®cos (3ni5 — 5a> + 252) 


^■( 114 ) 


cos (■» + «) — 252) = (1 — cos (nt + (o — 252) + le® cos (Znt — <o — 252) 

— -Je® cos (lit — Z(o + 252) — ^y^cos (Znt + <o — 452) 

+ ■|y®co8 (nt + 3w ~ 452) + c cos 2 (nt — 52) — e cos 2(o>— 52) 



cos (v — 3w + 252) = (1 — ^ cos (nt — Z(o + 252) — ■Je® cos (ni + a) — 252) 
+ Ic® cos (Znt — 5(0 + 252) + -J-y® cos (nt + 3a> — 452) 

— ■J-y® cos (nt — 5(0+ 452) + e cos (2nt — 4tt) + 252) 

— c cos 2 ((0 — 52) + ^y* cos (nt — (o) — ^-y® cos 3 (nt — (o) 


L ; (116) 


cos (Zv~co— 252) = (1 — 9c®) cos (Znt —eo— 2,52) + •®^c® cos (5nt — 3« — 252) ^ 
+ -^6® cos (jii + ® — 252) + |y* cos (3n.i + (M — 452) 

— fy® cos (5nt — « — 452) + 3c cos (4ni — 2a> — 252) 

— 3e cos 2 (nt — 52) + |y® cos (nt — w) — |y® cos 3 (nt — ai) ^ 


m 


If we now substitute these values in equation (105), and at the same time 
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change v into nt in the terms of the fifth order in that equation, it will 
become 


^ — 1 — e* + e {1 — •|e“ + -3^e^}cos(7i.i — <a) 

+ e* {1 — cos 2 {nt ~ay) + {|e» - cos 3 (nt-oi) ^ 
+ cos 4 {nt — <w) + cos 5 {nt — w) ^ 


(118) 


In the substitution in equation (105), the coeficienta of all the terms contain 
ing SI become identically equal to nothing, and the quantity y disappears from •* 
the coefficients of all the remaining terms, whence it follows that the radius 
vector of the orbit is entirely independent of the position of the orbit We shall 
therefore have 


+ cos {nt — w)+^{l— cos 2 {nt - a>) 

+ {I®* ~ cos 3 {nt — <y) + |e‘ cos 4 {nt — tu) + f|4e®cos 5 {nt — w) 


(119) 


In order to find the tangent of the latitude it is only necessary to find 
sin {v - Q,), and multiply it by y To find sin {v - Q) we shall observe that if we 
change cos {mnt + a) into ain(mnt + a), and —sin{mnt+a) into +cos(wni + a), 
cos (’^ + «) ■will change sin {mv + a) If we then make these changes in equa- 
tion (113) and put m = l, a = -Sl, we shall get, after multiplying by y, the 
following value of y sin {v - Si), or tan 6, 


tan<? =y {1 - e* — |.y* -{- ^ -p 

■“ {if ~ ilry® — liey } sin 3 {nt - ^3) + sin 5 {nt - Si) 

-fey {1 — fe® — ■J-y®}sin(2nt — — g3) — ey sin{a» — Q) 

+ {iff - fieV + -hf - sm (3?i< -2o)— Si) 

+ {if' -if - i^*y + •^y*+ ■^ey}sin {nt - 2cu -f £3) 

+ {if—hf-Hf^} sin {nt + 2oi- 3£3) sin 3 (u» - £3) 

-f fe®y sin (4»i — 3a) — £3) -f sin {4^t—Zo) — £3) 

+ { -W} sin {2m,t - 3a) -f £3) -f sin (5»«-4a) - £3) 

{ihsf' — i^f^} Sin {Bnt ~ 4a) -f £3) -f -j-lyy® sin {nt -f 4a) — 5£3) 

'I' {fif^ — xisf} sin {nt — 4a» + 3£3) — -ley* sin {4nt — <o — 3£3) 

-Ay*sin(3n.t-f2a)-5£3)-fiey*sin(2?i«-fa)-3£3)-|ieVsin(6nt-2a)-3£3) 


( 120 ) 



DEVELOPMENT OP THE ELLIPTICAL MOTION. 


31 


This equation will give 


J tan» ^ { 1 - 6“ ■ - sin {nt - S2) - { 1 — 9e® + |y } sin 3 (nt - S3) I 


+ sin 5 (nt — S3) — ^ey^ sin («) — S3) + sin (2nt + a> — 3S3) 

+ sin (2nt — <o — S3 )— ^ef sin (4:nt — a» — 3S3) - sin(37!,t +2a)-5S3) 

+ shi (3ni — 2£o — S3) — sin (6nt — 2a) — 3S3) 

+ 'sV?^ ~ /} sin (nt — 2 ( 1 ) + S3) + {/ — 7e®} sin (nt + 2a> — 3S3) 


^.( 121 ) 


•J-tan®0 = -|y*8in(nt — S3)— •]^y'sin3(nt — S3) H--^y®sin5(?it — S3). (122) 


Now we have ^ = tan 5 — J- tan* } tan ' d —, etc. (123) 

If we now substitute these values of tan^ and its powers in equation (123) we 
shall get the following expression for the latitude of the moon, which is correct to 
terms of the fifth order : 


^ = y {1 - e* - fy* + + ffeV} sin (nt — S3) 

“ {-sV/ ■“ tSt/ “ fK}'®} sin Z(nt — S3) + sin 5 (nt — S3) 

+ ey { 1 — I e® — ly*} sin (2nt — «» — S3) — ey { 1 — ^y*} sin (w — S3) 

+ {fc'y - f J6'V + — flfl*y*} sin (3n« - 2u» — S3) — •Jey* sin 3 (o) - S3) 

+ {i:«V - -If - sin (nt - 2w + S3) 

+ {iy® “ ‘isf ~ sin (nt + 2a) — 3S3) + |e’y sin (int — 3<u — S3) 

— ■Jcy’ sin (int — a* — 3S3t) + sin (2n.t + a» — 3S3) 

+ “ i®y®} sin (2nt — 3«> + S3) + l-lfeV sin (5nt — 4w — S3) 

■i" {x'i's^V “ •A'fi'y®} sin (3nt — 4w + S3) + sin (nt +■ 4a> — 6S3) 

+ {uVeV “ ri?/} sin (nt — 4« + 3S3) — -a^^y' sin (3?it + 2«> — 5S3) 

— ^fi*y®8in (6n,t — 2w — 3S3) 


. (124) 


7. We have thus found the values of the three co-ordinates r, v and ff in terms 
of the time. We may however find (9 directly from the differential equation (28), 
which will serve .as a verification of all the developments which have thus far 
been made in the determination of t, v, r and d. For this purpose we shall sub- 
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stitute tlie values of e, e' and o" in equation (28), by "wbicb means it will 
become 

Zcos(v-Q) (125) 


Now equation (119) will give 

■^ = 1 + + |e* + {2e + |e® + ffe®} cos (nt - o) 

/ 

+ {fe* + -^e^} cos 2 (nt — a>) + cos 4 (nt — w) 

COB 3 (n6 — eo)-i-^^>^e^ cos 5 (nt — a>) 


(126) 


and equation (113) gives, by putting m = 1, and a = — Q, 


cob(v — Q) = + coa(nt—Q) 

— e cos («) — ^3) + {e — + ^ef} cos {2nt — a) — Si) 

+ {K ~ Hfi* + Ay* + A«V} cos (3nf — 2a> — £S) 

— {if -Thy* -lieV} cos 3 {nt - Si) 

+ “Ay* ~ cos {nt + 2ai — 3£S) 

— ii^ + iy* - Ac* - Ay* ~ A A } cos {nt — 2co + Si) 

+ f e® cos {4M — 3(0 — Si) — {A®^ cos (2nt — 3w + £3) 

—^ey^ cos 3 (o) — £S) + cos {2nt + o — 3£2) 

—•fey* cos (4?^^ — co— 3Si) + cos {6nt — 4tt) — Q) 

~ {tA®* + A®V*} cos {3nt — 4a> + £3) + -jhf cos 6 (nt — Si) 

+ xAy* cos {nt + 4 : 0 ) — 5Si) — {jhf + A®V} cos {nt — 4<u + 3£3) 
— Ay* COB {3nt + 2£i> — 5£^) — flA* cos {5nt — 2co — 3Si) 


(127) 


1 

If we now multiply equations (126) and (127) together, and the product by y, 
we shall obtain 
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= 1 {y-^e^y + 1/ -*eV - AeV}cos (nt - Q) 

r cr 

~ {iy* ~ ll^y* — M®V} cos 3 {nt — Si) + jf^y* cos 5 {vi — Si) 

+ {2ey — + \ef] cos {2nt — <d — S3) + -^eV cc® — 3<u — S3) 

+ + 'H'CV + ift'/} cc® — 2a> — S3) 

+ {^y® ~ -hf “ iicy } CO® + 2ft> - 3S3) 

+ {ie*y “ if + ^cV + ^y' + cos (nt — 2w + S3) 

— {^ey^ — cos {2nt — 3a» + S3) + cos {2nt + 6) — 3S3) 

~ McV} cos {Snt — 4a) + S3) — ^ef cos (4ni — o) — 3 S3) 

■*■ {'b'jcV "liry*} cos — 4a) + 3S3) + cos (6ni — 4a) — S3) 

+ xf-g-y* cos {nt + 4a) — 6S3) — cos (3n4 + 2a) — 5S3) 

— -f-leV cos {5nt — 2a) — 3S3) 

If we multiply this equation by 

^/af^VlzI = T/a^{i_^e?-iZ-ieH|Z + i6V, 

/l + Z 

and put we shall obtain 

y [I- l-Z + + HcV} cos {nt - Si) 

ndt 

~ ^y^ { ^ “ My* ^®'C*} cos 3 {nt — S3) + xfvy* cos 5 {nt — S3) 

+ 2e)'{l — |e*“|Z} cos(2?it— 0)— S3) 4-'Y6*ycos(4ni — 3o) — S3) 

+ {'V®*y “ “ f 1<5V + T^jZ) cos {Bnt — 2a) — S3) 

+ {^y*_.j;?jj,»-^|ey}cos(n«+2^«)--3S3)-||cycos(5»i-2a)-3S3) I (129) 

{i^ ~ ic* + TsVy* ■^jc*y*} cos (wt — 2f«) + S3) 

— M-eZ “ cos {2nt — 3(o + S3) + ^-eZ cos {2nt + o) — 3S3) 

4. — fleV } cos {3nt — 4w + S3) - -J-eZ cos {int — co — 3S3) 

+ {-B^cV - rbZ} cos {nt - 4w + 3S3) + cos {5nt -la — Si) 

+ xfirZ cos (wt + 4:0)— 5S3) ~ •g*^y* cos {3nt + 2w — 6S3) 
r 
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Equation (129) gives, by integration, 

d = {6) + -^e* + ^/ + f|eV}sm(n,i-£3) 

- {^7® - rhf - fi«V} sin ^ {nt - ^) + sin 6 (nt - Q) 

+ 6y {1 — sin (2nt — ai — ^) — ^eVsin (5nt — 2(o — 3SJ) 

+ {|eV - HeV + */ - W} sin (3nt -2a>-Q) 

+ - iy* “ -h^^y + -hf + sin {nt — 2o}-{-^) 

(i/ ■" ~ W®V*} sin -\-2ci) — 3Q) + -feV sin (4ra5 — Sco — Q) ' 

— -^ey* sill (4nt — ai — 3g3) + ^ey® sin (2nt + <u — 3gS) 

{t^®V “ i®y*} sin (27ii — 3ai + sm — 4tt) — £2) 

+ {rl^sV ~ sin (3n« — 4ai + £S) + xls-y® sm (nt + 4ai — 5£2) 
b {■^^V ~ tIt/} sm — 4a) + 3SS) — -^y® sm (3wt + 2a) — 5£2) 


(130) 


(d) being the constant quantity to complete the integral It is evident that the 
two expressions for the value of d given in equations (124) and (130) will be iden 
tical if we make the constant 

(d) = — ey {1 —\f} sm (w — Q) —-ley® sin3 (w — Q, (131) 


■which makes equation (130) satisfy the condition that the latitude of the moon 
shall be equal to the latitude of the perigee of the orbit when nt = co, or when 
the moon is at the extremities of the transverse axis of the oibit The perfect 
agreement of these two determinations of the value of d proves conclusively that 
all the preceding analytical developments have been correctly made 

8 In order to show, by a few numerical examples, that the values of v, r and 
6, given by equations (104), (119) and (124), are correct, and at the same time 
show that the development of these quantities m senes to terms of the fifth order 
IS sufficient, m the lunar theory, 5je shall now reduce these equations to numbers 
by using the values of e and y, corresponding to the elements of the moon s orbit 
which were employed by Delaunay in reducing his equations of the lunar theory 
to numbers These values are, e = 0 06489930, and y = 0 09004660 And 
instead of giving the value of r in terms of the moon’s mean distance as the unit, 
we have multiplied equation (119) by the constant term of the moon’s parallax, 
supposing it to be equal to 3422 "3 
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We shall therefore find 

V = nt + 22638."97 sin {nt — co) + 777."074: sin 2 (jit — m) 

+ 36. "997 sin 3 (yit —lo) + 2."010 sin 4 (nt—(o) -f 0."118 sin 6 (nt—w) 

— 411."374 sin 2 [nt - gj) + 415."469 sin 2 (w - £3) 

— 45."255 sin (3nt — (o — 2£3) + 45."691 sin {nt + co — 2Q) 

— 4."096sin (4nt — 2W — 2g3) +0."093 sin(wt + 3a> — 4£3) 

— 0."186 sin (Znt + <o — 4g3) + 0."093 sin (6?it — (o — 4£3) 

— 0."340 sin {5nt — 3co — 2g3) — 0."848 sin (%it + 2co — 4g3) 

4- 0."424 sin 4 (nt - £3) + 0."424 sin 4 («» - £3) 

d = 18461."23 sin {nt - £3) + 1012."716 sin {2nt - a» - £3) 

- 1017."590 sin {to - £3) + 62."619 sin {Znt - 2<o - £3) 

+ 18."584 sin (nt + 2w — 3£3) — ll."662 sin {nt — 2oj + £3) 

-- 6."993 sin 3 {nt — £3) + 2."067 sin {2nt + co — 3£3) 

- 1."0335 sin (4»t - w - 3£3) - 1."0335 sin 3 {co - £3) 

+• 4."098 sin (4nt - 3<<^ - £3) + 0."2746 sin {5nt - 4^; - £3) f' 

+ 0."0286 sin {nt + 4w — 5,Q) — 0."1206 sin {5nt — 2co — S£3) 

— 0."0191 sin {Znt -I- 2co — 5£3) ~ 0."0024 sin {nt — 4«; + 3£3) 

+ 0."0057 sin 5 {nt - £3) - 0."7773 sin {2nt - Zco + £3) 

— 0."0520 sin {Z?it — 4« + £3) > 

== 3422."300 + 187."811cos(n<-w; +10."304 cos2(nt-fi>) ' 

+ 0."636 sin 3 {nt — co) + 0."0414 cos 4 {nt — co) >■ ; 

+ 0."0028 cos 5{nt — co) 

7T denoting the moon’s parallax. 
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We shall now compare these approximate values of v, d and n with the values 
derived by an exact calculation from rigorous formulae for the same co ordinates 
For this purpose we shall observe that if we denote by u the eccentric anomaly, 
and by nt — co' the mean anomaly measured on the plane of the orbit, the rela 
tion between u and nb — co' will be given by the equation 


nt — 0 )' = u — esmu 


(135) 


Then we shall have the following equation for the deternunation of the true 
anomaly v' — co', and radius vector r. 


tan ^(-b' — (o') = ^ tan 


l/l-e 


r — a(l — e cos v) 


(136) 


Equations (135) and (136) are entirely rigorous , and if, for a given value of 
the mean anomaly nb — co', we compute the corresponding value of u with accu- 
racy by means of equation (135), we can then compute the corresponding values 
of u' — co' and r by means of equations (136) Now we also have the value of 
a>' — SS by means of the equation 

tan {co' — Si) = tan {co — Q) sec (137) 

then adding the value of v' — co' to co' — we shall have the true distance 
v' — Si, of the moon from the node measured on the orbit 

Then the spherical triangle vv'Si (see figure. Art 10, Introduction) will give 


and 


tan {v — Si) = tan {v' — Si) cos i 
sin d = sin (u' — £2) sin % 


(138) 


The formulae (135-138) will give the rigorous values of r, v and d correspond- 
ing to any value of the mean anomaly nt — co^ 

If we suppose that 5° 8' 43 "28, and — = we shall find that 
o)' — Q, = 30° 6' 0 "99 Then if we designate the values of v and d, derived from, 
equation (138), as "observed values, and the values of the same quantities 
derived fiom equations (132) and (133) as "calculated values,*' we shall obtain 
the following table showing the values of the co oidinates v and d for the values 
of the mean anomaly given in the first column of the table 
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Mean 

Anomaly. 

V — 

Obsei 

■Q, 

WEE. 

6 

V — 

Calcu 

Q 

— 

LATBD. 

e 

Calculated.- 

Av 

—Observed. 

Ad 

30° 

o 

63 

20 

58"23 

o 

+ 4 

36 

4(^30 

63° 

20' 

5^^06 

o 

-i-4 

36 

430 

-d'i7 


40 

V4 

18 

13.61 

4 

57 

16.05 

74 

18 

13.72 

4 

57 

16.06 

-1-0.11 

4-0.01 

60 

95 

45 

20.32 

5 

7 

10.39 

95 

45 

20.51 

5 

7 

10.34 

-fO.19 

-0.05 

80 

116 

27 

0.21 

4 

36 

33.14 

116 

27 

0.22 

4 

36 

33.10 

+ 0,01 

-0.04 

100 

136 

19 

36.09 

3 

33 

29.25 

136 

19 

35.97 

3 

33 

29.18 

-0.12 

-0.07 

130 

164 

46 

9.63 

-1-1 

21 

18.38 

164 

46 

9.80 

-l* 1 

21 

18.40 

4-0.17 

4-0.02 

160 

192 

4 

24.48 

-1 

4 

44.42 

192 

4 

24.44 

-1 

4 

44.40 

4-0.04 

4-0.02 


The columns of residuals show that the development of the longitude and lati- 
tude in series as far as terms of the fifth order can in no case be in error to a 
greater extent than about 0."2 in longitude, and 0."1 in latitude; and it would 
therefore seem that the development of the lunar theory to that degree of approxi- 
mation would be amply sufficient for all the purposes of astronomy, except for 
those terms of perturbations producing the inequalities of long period. We would 
here also observe that the value of the parallax given by equation (ISd) can in no 
case differ from the rigorous value of that co-ordinate to a greater extent' than 
O/'l ; and a more perfect agreement between approximative and rigorous formulae 
could hardly be expected or desired. 

9 . The whole of the mathematical theory of the moon’s motion is contained in 
the preceding articles (1-8), when we neglect the consideration of the effects 
of the disturbing forces. The expressions for the three co-ordinates v, d and r, 
given in equations (104), (124) and (119), are perfectly homogeneous in their 
development, every circular function or argument of the different terms of the 
equations being measured on the fixed plane of projection instead of the plane of 
the orbit. There are therefore no terms depending on the mean distance of the 
moon from the perigee measured on the plane of the orbit, usually called the 
mean anomalies; but the equivalent terms depending on the difference of the 
mean longitudes of the moon and of the perigee measured on the plane of projec- 
tion is given instead. For this reason it is unnecessary to make any reference to 
the distance of the perigee from the node of the orbit measured on the plane of 
the orbit, but simply to designate the longitude of the perigee by the distance of 
its projection from the origin of longitudes and measured on the fixed plane. And 
therefore the expressions for these three co-ordinates which we have given are 
simpler than the equivalent expressions given by Delaunay, who reckons longi- 
tudes on the fixed plane,- and the other arguments of his equations upon the plane 
of the orbit, referring them to the fixed plane in terms of certain angular func- 
tions depending on the distance between the node and perigee, and also to the 
mean anomalies which are measured on the plane of the orbit, Delaunay’s 
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theory therefore requires the consideration of two planes of reference, one for the 
longitudes and the other for the anomalies For this reason his formulae are less con- 
veniently adapted for developing the perturbations than the expressions which we 
have given We may observe, however, that Delaunay s expressions for v and d 
may be made identical with equations (104) and (124), as far as terms of the fouith 
order, by changing the y of his notation into •i'7(l“~f7^), and also changing g into 
(w — SS) + i/^Bin 2(a) — Q>) We may then correctly put his l^nt — w, and then 
the two systems of equations will be made identical Considering the intiicacy 
and importance of the subject, and the great extent to which it is nf^ccssary to 
cairy the approximations, it is very essential to reduce the fundamental elements 
of the problem to as simple a form as possible , and it is believed that the differ- 
ential equations of the moon^s co ordinates, v, 6 and r, given in § 1, leave little to 
be desired so far as the determination of these quantities is concerned We ^^hall, 
however, in order to render this chapter more complete as a basis for the further 
development of the lunar theory, now give the analytical expressions for the vari- 
ations of the elements y, fj, a, e and a>, depending on the action of the disturbing 
foices 

10 For this purpose we shall observe that we have already found 

tan ^ (139) 

If we suppose thab o' and c" are variable, from tbe action of the disturbing 
forces we shall get, by differentiating equation (139) and dividing by dt, 


dQ> _ cos^Q ^de'' dG' '\ 

dt d t dt d dt J 


(140) 


But we have c' = cy cos £3, and 
become 


dt ory 


e" = cy sin ^ , therefore equation (140) will 


cZc" Bin Q do' 


dt 


cry dt 


(141) 


If d and c" are variable by reason of the disturbing forces, equations (9) and 
(10) will give 
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If we now substitute these values of and in equation (141), and also 

at at 


put 


we shall find 


0 = 


— (1 < 


yi + f 


vT+7 


dQ 


dt y afi (1 


^ (, 44 ) 


In the same manner we may compute the differential variations of the 
elements y, e and coj but the process is simple, and it is unnecessary to give 
it in detail We shall therelore insert here the equations to be employed 
in the further development of the lunar theory, the solution of which cannot 
fail to give the moon s co ordinates with all the precision required by obser- 
vation, if the inequalities of her motion are produced by the forces of gravi- 
tation 

If we therefore substitute the values of c, c', c", / and in equations (19), 
(20) and (21), we shall obtain the following differential expressions of the co ordi- 
nates Vj d and r 

in. = J 1 _ + / 

dt i/l r*cos*5 (. — ^ J\dv/ )’ 


dt 


( 0 ) 
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^ ^ecos^dcoa^sinfa — o’)!! Tf^V^l 

dt VafiiX^T 1 i/a/.(l-e?)‘^\^W i 

_ cos 6 cosvj^ |r*{cZ^ mid smv — dv cos0 coav} 


+ 


l/a/i(l-e“) 


+ { 2 r(i'ucoa^smv ■ 


. 55 £i!cZr}/^ 

cos^ \dv J 


■ — cos d sm V 

l/afiil-e^ 


+ {2rd0 coa dBinv + dr sin 6 sin v} 

)J |r*{cZ^ sm d cos-y + dv cos d sinv} 
a-y , ^ [(. 


+ {2rdv cos d coBV ■ 


/^\ 

ydv] 


+ {irdd cos <? COS'D + dr sin 6 co&v} | 


t/T+?' 


'l/afi (1 — 


S{T)- 


(D) 


dt -i/afiil-^l \d'vJ y \<^^/) 




dB\^ 


fiecoadt^^ =r® |sin^sin(i? — <u)^ — cos5cos('D — <y)^ (dr/ 

, fcos('D— (w) dr Of -.d'u) /di?' 

H-J — i — — 2rcosDSin('D — w)— W — 

1 cos^ dt dt)\dv 

f ad'>' <n / v/dJ2' 

—4 sin^— + 2r cos D— y sin (v — «) — - 
1 dt dij \d^; 


(B) 


(B) 


(G) 
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— cos/?osin<? cos(^)-w) }^ — cos^o cosd sin(i; — (~j 


{ 


cos/9« 


.dr 


siu {v — + 2r {cos cos d cos (« — <y) + sin sin ^ } — 1 Z' 

cosd? dd ^ 


+ 


dr 

dt 


L(H) 


dR 


{sin 6a cos 6 — cos 0^ sin 0 cos (u — O))} 

— 2r {cos^oC 08 ^coa('!; — <u) + sin^osin^}— f 

dt 


And lastly, equation (18) will give 



Wo have already integrated equations (B), (0) and (D), when the function iJ 
is equal to nothing; and if we substitute these approximate integrals in equations 
(E-I), wo shall obtain the approximate values of the variations of the elements 
of the orbit. Wo must then substitute the corrected elements a, e, y, co and Q, 
in equations (B), (C) and (D), by which means we shall obtain corrected values of 
V, 6 and r, which will contain terms depending on the first power of the disturbing 
force; then with the corrected values of the co-ordinates v, 0 and r, together with 
the now elements of the orlht, wo must repeat the solution of equations (E-I), by 
which moans wo sliall obtain corrected values of the variations of the elements, 
with wliich wo can repeat the solution of the equations which determine the 
co-ordinates v, 0 and r. By continuing the solutions in this maimer wo may 
obtain all the terms of sensible magnitude arising from the disturbing force. 

11 . In bringing to a claso this first chapter of a lunar theory, it will be inter- 
esting to compare the results of our amalysis with the conclusions arrived at by 
other mathematicians, in so far as the different methods of investigation will 
permit. For this purpose we shall ob.sorvo that both La Place and Plana in 
their lunar theories have first determined the mean longitude of the moon in 
functions of the true longitudo as tlio independent variable, in the same manner 
as wo have done in §4 of this chapter; and as the co-ordinates of the funda- 
mental differential equations employed by them have precisely the same signifi- 
cance as in this olnapter, it would seem that the several results arrived at should 
be. idcnticiilly the same, provided the methods of development employed were 
equally goiit'rnl and the calculations had been correctly made. 

La Pi.acr gives the two following equations for the determination of the 
moan longitude and tho reciprocal of the radius vector, in functions of the 
8 
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true longitude as the independent variable, whicb are independent of tlie 
perturbations 

wt = u - 2e (1 - if) sin (u - to) + |e^ sm 2 (u - to) ~ sin 3 (u - to) 

+ if sin 2 (r - gj) - ^ Bin (3u - to - 223) - f ef sin (u + to - 2g5) 

-=l + e® + e(l + e®-if)co8(u-to) + iey® cos (3u - to - 2gJ) 
r 

+ ief cos (u + to — 2Q) 

■while equations (89) and (105) mil give, to terms of the same order, by neglecting 
the constant introduced by integration. 




rA = u — 2e sin (u — to) + f e® sin 2 (u — to) — ie® sin 3 (u — to) 

+ if sin 2 (« — g^) — ief sin (3u — to — 2gS) — if sin (u - 3to + 2gS) 


(147) 


- = l + 6® + 6(l + e®)cos(u — to) +ief cos(3d — to — 2g5) 
r 

— iey® cos('U + to — 2g3) + iey cos (‘0 3to + 2g3) 

Comparing now the values of nt given by equations (145) and (147), we notice 
that the coefficients of sin (u - to) in the two equations differ by the quantity Jef , 
or by a term of the third order We also see that La Place s equation contains 
a term of the third order depending on sin(u + to — 2Q), which does not appear 
in our equation , while our equation contains a term ot the third order depending 
on sin('o-3to + 2g2), which does not appear in the equations of La Place or 
Plaka except with a coefficient of the fifth order These two last mentioned 
differences arise directly from the introduction of the latitude of the perigee into 
our analysis, an element that has never been explicitly considered m the analysis 
of the lunar and planetary theories, but which will be found to affect all teims of 
the third and higher orders in the planetary theories to the same extent as in the 
lunar theory Indeed, the usual form of development is such that the aigument 
of an inequality indicates the analytical order of magnitude of its coefficient, 
because the quantity a) is always connected with e, 2io with e®, etc , while g^ is 
always connected with y, 2gS with f , etc , so that the combination of the terms 
containing the sines and cosines of angles into which these quantities enter, pro- 
duce coefficients with powers and products of e and y in which the sum of the 
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indices is equal to tlie arithmetical sum of the coefficients of o> and Q in the 
expression of the angle But the introduction of the latitude da, of the perigee, 
brings o) into the analysis unaccompanied by e, and in such manner that we have 
cosdo = {1 + y*sin®(<w — which, being developed, produces the second 
member of equation (77), an expression in which the order of the coefficient of 
any cosine is only one half as great as the sum of the coefficients of w and Q, in 
the given angle 

It IS therefore easy to see that the product of two or more sines or cosines 
having coefficients whose order of magnitude has different relations to the con 
stant terms of the angles will completely vitiate the law in regard to the order of 
magnitude of any coefficient which would take place, provided the coefficients of 
all the sines or cosines bore the same relation to the constant terms of the angles 
And hence the term depending on sin {v + m — 2£S), of La Place s and Plana s 
analysis, disappears, and the term depending on sin (u — 3<w + 2Si) comes in with 
a coefficient of the third order 

If we now suppose that ai = gj, in equations (145) and (147) they both reduce 
to the following 


v = nt — {2e + sin (v — w) + sin 2 (■» — oi) — (-le® + ^ef) sin S(v — to) ) „ 

> (i4y) 

+ -J/sin2(u — g^) J 

In comparing equations (146) and (148), we notice that the coefficients 
of coa(v — <o) differ by the term ief, while equation (146) has the term 
^ey®cos(u + — 2£S), instead of the two terms, — Jey® cos (wt + tu — 2g5) 

+ cos (v — 3<o + 2Q), in equation (148) However, if we suppose that = g3, 
they both reduce to the following 

^ = 1 + e® + e (1 + e® — Jy®) cos (v — w) -j- Jey® cos 3 (u — cc) (150) 

If we suppose, however, that v = o>, in equations (146) and (148), they will 
become respectively 

^ == 1 + e + e® + e® — + Jey®cos2((w — gS), (151) 

and - = l + e + e® + ^ (152) 


a 

r 


But the true elliptical value of -, when the moon is in perigee, is 

CL ^ 

■ — = 1 + e + e® + e®, which is the same as equation (152) , and therefore 
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equation (151) is erroneous, unless o) = Q, in ■wliicli case it becomes tke same as 
equation (152) If we now take tbe differential coefficient of equation (146), we 
shall obtain 


=e(H- ^ - iy*) sin(u - ci»)+ f ef sin(3u - oi - 2g5) + Irey® am(v+eo-2Q,) (1 53) 

T^dv 


J 

It is evident that this value of — never becomes equal to nothing when v = co, 

dv 

unless (0 is also equal to £3, or to £3 + 90°, a conclusion already deduced from 
equation (16) of the Introduction This defect, however, does not exist in equa 

J 

tion (146), for it gives — =0, whatever be the relation between 

dv 

the quantities co and Q, 



